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FOREWORD 


The  ARIES  model  was  developed  by  Calspan  Corporation  (formerly 
Cornell  Aeronautical  Laboratory,  Inc.)  as  one  task  of  its  System  Engineering 
Technical  Assistance  (SETA)  contract  with  the  AGM-86A  Program  Office  (RW  86) 
of  the  Aeronautical  Systems  Division  (ASD)  of  the  Air  Force  Systems  Command 
(AFSC).  The  model  was  defined  during  the  period  February  1972  to  Juie  1972 
wder  Contract  No.  F33657-72-C-0228.  Portions  of  the  model  were  programmed 
on  the  IBM  370/16S  conputer,  located  at  Calspan,  Buffalo,  New  York,  during 
the  period  July  1972  to  Septeafcer  1972  under  Contract  No.  F33657-72-C-1013. 

The  final  testing  and  checkout  of  the  model  computer  program  were  accomplished 
on  the  CDC  6600  computer  at  ASD,  Wright-Patterson  AFB,  Ohio. 

This  report  contains  a general  discussion  of  the  scope  and  practical 
applications  of  the  model  from  the  Integrated  Logistics  Support  (ILS) 
manager's  point  of  view,  a rigorous  technical  discussion  of  the  model  for 
operations  researchers,  and  a description  of  the  model  computer  program  and  its 
applications. 

Lt.  Col.  L.  Davila- Aponte  and  Mr.  John  Burchett  of  the  ILS  Division 
of  the  AGM-86A  Program  Office  were  the  USAF  monitors  during  the  development 
of  the  model. 

Messrs.  E.  Pringle  and  T.  Wojcinski  and  Dt.  W.  Fryer  contributed  to 
the  model  computer  program  development  at  Calspan' s Buffalo,  New  York, 

Advanced  Technology  Center,  and  Dr.  Sol  Kaufman  of  Calspan 's  Operations 
Research  Department  at  the  Center  provided  technical  assistance  in  the  review 
and  development  of  the  analytic  structure  of  the  model. 


ABSTRACT 


V-_— 

managers  in  their  planning.  It  enables  then  to  optimize  recoverable  item 
conponent  reliabilities,  maintenance  capabilities,  and  maintenance  strategies  - 
as  well  as  initial  stockage  levels  for  the  multi-echelon  base-depot  supply 
system  --  while  taking  into  account  the  corresponding  impact  upon  system  avail- 
ability and  life  cycle  cost  for  unscheduled  maintenance.  > 


^Item  demand  is  assumed  Poisson.  There  are  four  alternative  main- 
tenance strategies  for  an  item,  and  the  strategy  choice  is  a variable  in  the 
problem.  Rather  than  considering  only  an  item's  unit  cost,  the  model  confutes 
the  item's  10-year  life  cycle  cost  f (^unscheduled  maintenance.  The  sum  of 
item  expected  backorders  across  aU/oases  is  minimized,  subject  to  a life 
cycle  cost  constraint  for  unscheduled  maintenance.  The  item  expected  backorder 
versus  cost  function  can  lvatfe  large  regions  of  non- convexity.  The  resource 
expenditure  algorithm  for-  the  model  incorporates  allocations  from  these  non- 
convex  regions  te  provide  solutions  near  a constraint.  The  maximum  region  of 
uncertainty^for  the  minimized  expected  backorder  sun  is  defined,  and  this 
provide*  a basis  for  judging  the  quality  of  the  result. 


ARIES  evaluates  the  minimized  expected  backorder  solution  to  deter- 
mine the  expected  number  of  NORS  vehicles  for  both  the  case  where  conqplete 
cannibalization  of  parts  is  assumed  and  the  case  where  there  is  no  cannibali- 
zation. The  ARIES  solution  defines  the  expected  backorder  value  for  each 
item;  this  parameter  is  a necessary  input  to  the  MSCAM  model  to  perform 
refined  reliability/maintainability  trade-offs  on  the  individual  items.  The 
manpower  and  item  stockage  level  requirements  for  the  individual  bases  and 
depot  are  determined  by  the  ARIES  solution. 
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A weapon  system  designed  to  perform  near  the  state-of-the-art  in  a 
performance  sense  (e.g.,  in  range,  speed,  accuracy)  is  a dubious  asset  if  it 
is  nearly  always  down  because  of:  (1)  equipment  failure,  (2)  a low  probability 
of  repair  within  an  acceptable  time  interval,  or  (3)  a high  probability  of  a 
shortage  of  the  replacement  spare  parts.  Consequently,  "availability”  becomes 
an  important  measure  of  total  system  effectiveness.  Availability  is  the 
resultant  interaction  of  the  quantitative  aspects  of  system  component  reliabil- 
ities, maintainability  capabilities,  and  supply  support  effectiveness. 

Additionally,  it  is  a design  goal  to  achieve  an  acceptable  level  of 
system  availability  subject  to  minimizing  the  total  life  cycle  cost  for  the 
weapon  system.  The  total  life  cycle  cost  for  an  item,  segment,  or  system  is 
the  sum  of  acquisition,  operation,  and  maintenance  costs.  This  report 
discusses  a mathematical  model  that  assists  the  PO  in  the  decision-making 
process  by  deriving  system  conponent  reliabilities,  maintenance  strategies, 
and  deployment  stockage  levels  while  taking  into  account  the  corresponding 
impact  upon  system  availability  and  life  cycle  costs  for  unscheduled  main- 
tenance. 

The  technical  discussion  of  the  ARIES  model  (Section  IV)  was  developed 
and  written  in  a style  suitable  for  submittal  to  a professional  journal  for 
review  and  potential  publication.  The  journal  requires  that  papers  be  double- 
spaced, reference  rather  repeat  previous  work,  etc.;  therefore,  to  save  the 
time  and  costs  associated  with  re-doing  the  technical  discussion  of  the  model 
in  the  typical  Calspan  report  format.  Section  IV  is  essentially  an  unaltered 
version  of  the  "journal-style"  paper.  The  double- spacing  offers  the  reader 
additional  freedom  to  make  notes  and  comments  within  the  text  of  the  dis- 
cussion. 
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The  remaining  sections  of  the  report  were  prepared  using  the  con- 
ventional reduced  line  spacing  but  in  a format  that  is  consistent  with  Section 
IV  to  assist  in  the  ease  of  reader  comprehension.  The  General  Discussion  in 
Section  III  has  been  developed  with  the  objective  of  minimizing  the  amount  of 
technical  details  while  still  providing  the  reader  with  a sufficient  overview 
of  the  model. 


Section  21 
SUMMARY 


This  report  presents  a technical  discussion  of  a conputerized  mathe- 
matical model  developed  by  Calspan  Corporation  to  assist  the  AGM-86A  Program 
Office  in  evaluating  and  in  reducing  the  logistics  costs  of  the  SCAD  weapon 
system.  The  model  is  applicable  to  other  weapon  systems  that  have  a logistics 
system  which  is  similar  to  the  SCAD  base-depot  concept. 

The  ARIES  model  is  used  to  derive  the  optimum  weapon  system  logistics 
configuration  for  unscheduled  maintenance  for  recoverable  items  (Shop  Replace- 
able Units,  SRUs)  as  a function  of  system  availability  and  10-year  life  cycle 
costs. 

The  features  of  the  model  that  assist  the  ILS  manager  in  identifying 
and  reducing  weapon  system  life  cycle  costs  for  logistics  are: 

e Consideration  of  four  alternative  maintenance  strategies; 

e Calculation  of  life  cycle  costs  (unscheduled  maintenance); 

e Creation  of  a permanent  data  base; 

e Permitting  of  analyses  at  various  cost  constraints  without 
recalculating  the  data  base; 

e Use  of  simple  Poisson  distribution  fbr  demand; 

e Provision  of  a resource  allocation  algorithm  that  derives 
solutions  which  are  within  1 percent  accuracy  for  most 
realistic  problems; 

e Relation  to  operational  measures  of  effectiveness; 

• Provision  of  continuity  with  the  individual  SRU  reliability/ 
maintainability  trade-off  data  and  model; 

e Computation  of  base  and  depot  unscheduled  maintenance 
manpower  requirements;  and 

e Provision  of  useful  base  simnury  and  item-depot  summary  output. 
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Section  IQ 

GENERAL  DISCUSSION 


This  section  presents  a general  discussion  of  the  requirement  for 
and  the  practical  appli cation  of  the  ARIES  model  from  the  ILS  manager's  point 
of  view.  The  mathematical  details  of  the  ARIES  model  are  discussed  in  Section 
IV- 

ILS  MODELING  OBJECTIVES 

The  military  services  have  been  under  increasing  pressure  in  recent 
years  to  estimate  and  reduce  the  cost  of  developing,  deploying,  and  operating 
new  sophisticated  weapon  systems.  The  "cost  growth"  phenomenon  of  past 
weapon  systems  has  brought  about  enough  attention  that,  once  the  requirements 
have  been  established  for  a new  weapon  system,  the  following  order  of  priori- 
ties tends  to  prevail:  (1)  cost,  (2)  performance,  and  (3)  schedule. 

The  logistic  support  costs  for  a weapon  system  account  for  a signifi- 
cant portion  of  the  total  weapon  system  cost  over  the  life  cycle  of  the 
deployed  system.  In  1964,  the  Department  of  Defense  issued  Directive  4100.35 
describing  the  Integrated  Logistics  Support  (ILS)  concept.  Since  then,  the 
Services  have  been  required  to  consider,  estimate,  and  evaluate  the  life  cycle 
costs  (LCC)  implied  by  the  design  alternatives  encountered  throughout  the 
acquisition  process.  By  providing  improved  methods  for  making  trade-off 
evaluations  early  in  the  weapon  systems  design  phase,  attention  can  be  focused 
upon  reducing  life  cycle  costs  before  significant  funds  are  committed  to  a 
design  or  concept. 

It  has  been  the  objective  of  the  Calspan  technical  support  effort  for 
the  ILS  function  of  the  SCAD  program  to  review  previously  developed  ILS 
models,  to  modify /develop  appropriate  ILS  models  for  SCAD,  and  to  assist  in 
the  implementation  of  the  models.  The  MSCAM  model  (originally  developed  at  the 
RAND  Corporation  and  subsequently  modified  by  Calspan)  has  been  made  available 
to  the  SCAD  PO  to  evaluate  the  reliability /maintainability  characteristics  of 
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individual  SRUs  (Shop  Replaceable  Units)  and  how  these  characteristics  affect 
the  life  cycle  cost,  stockage  levels,  and  alternative  maintenance  strategies. 
The  ARIES  model  has  been  developed  by  Calspan  for  the  SCAD  PO  to  evaluate  the 
aggregate  effect  upon  system  availability  and  life  cycle  costs  for  all  SCAD 
SRUs. 

LIMITATIONS  WITH  PREVIOUSLY  DEVELOPED  ILS  MODELS 

A review  has  been  made,  and  is  continuing  to  be  made,  of  previously 
developed  ILS  models  that  have  been  suggested  for  application  to  the  SCAD 
program  by  organizations  other  than  Calspan.  The  following  paragraphs 
summarize  the  principal  limitations  of  some  of  the  previous  ILS  models  and 
enphasize  the  need  that  existed  for  a new  model  which  could  improve  the  SCAD 
ILS  modeling  applications  to  reduce  life  cycle  costs. 

ORLA  (Optimum  Repair  Level  Analysis)  Model 

ORLA  does  not  employ  a stockage  level  computation  procedure  that 
permits  the  planner  to  estimate  the  stockage  cost  required  to  meet  a specified 
level  of  system  perfoxmance.  Additionally,  the  effectiveness  of  computed 
stockage  levels  is  not  made  explicit.  Thus,  the  costs  computed  for  alternative 
repair  strategies  do  not  necessarily  correspond  to  the  same  levels  of  system 
performance.  The  ORLA  model  does  not  include  an  alternative  maintenance 
strategy  for  base-depot  repairs. 

AFLC  Drone- ISC  (Life  Support  Cost)  Model 

The  LSC  model  assumes  that  the  initial  stockage  level  is  equal  to 
the  peak  demand  month.  This  assumption  is  arbitrary  and  does  not  apply  to 
SCAD,  for  there  is  no  peak  demand  month  as  there  is  potentially  in  an  aircraft 
weapon  system.  The  model  does  not  consider  and  coapute  optimum  maintenance 
strategies. 


*The  terms  "SRU"  and  "item"  are  used  interchKgeably  throughout  this  report. 


METRIC  (Multi-Eche Ion  Technique  for  Recoverable  Item  Control)  Model 

The  METRIC  model,  developed  at  the  RAND  Corporation,  provided  the 
original  concept  for  evaluating  all  recoverable  items  in  a weapon  system  on 
an  aggregated  basis  while  minimizing  the  sum  of  total  expected  backorders 
across  all  bases.  However,  METRIC  has  the  following  limitations:  (1)  life 

cycle  costs  are  not  considered,  (2)  alternative  maintenance  strategies  that 
could  vary  with  the  level  of  LCC  expenditure  are  not  considered,  (3)  great 
emphasis  is  placed  upon  describing  the  demand  distribution  (i.e.,  Bayesian 
estimation  procedures  for  the  compound  Poisson  distribution,  which  entail  cal- 
culations that  are  of  questionable  marginal  value  considering  the  likely  accu- 
racy of  the  total  set  of  input  parameters),  (4)  the  item  condemnation  rates 
are  assumed  to  be  zero  (i.e.,  every  item  is  assumed  to  be  repairable  at  either 
the  base  or  depot),  and  (5)  the  output  data  are  not  as  readily  usable  and  com- 
plete as  is  desirable  for  the  SCAD  systems  analysis. 

GENERAL  DESCRIPTION  OF  ARIES 

ARIES  is  a computerized  mathematical  model  used  to  derive  the  optimum 
total  weapon  system  10-year  life  cycle  cost  for  unscheduled  maintenance  as  a 
fimction  of  system  availability.  The  operational  measure  of  system  availabil- 
ity is  the  expected  number  of  NORS  (Not  Operationally  Ready  - Supply)  missiles 
for  the  case  when  complete  cannibalization  of  parts  is  assumed  and  for  the  case 
when  cannibalization  is  not  permitted.  ARIES  is  applicable  to  any  weapon 
system  that  utilizes  the  base-depot  maintenance  concept. 

Within  ARIES,  the  measure  of  effectiveness  "the  sum  of  expected  item 
backorders  across  all  bases"  is  minimized  in  the  optimization  process.  A 
demand  for  an  item  that  cannot  be  satisfied  from  base  supply  is  a backorder 
or  due-out.  The  expected  backorder  measurement  takes  into  account  the  duration 
of  an  unsatisfied  demand  as  well  as  the  fact  that  it  occurred.  The  process 
of  minimizing  the  sum  of  expected  backorders  provides  a close  approximation 
to  the  process  of  minimizing  the  operationally  oriented  NORS  measure  of 
effectiveness  (no-cannibalization  case). 
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The  input  data  to  ARIES  include:  all  the  individual  SRU  data,  which 

are  essentially  the  same  data  processed  individually  with  the  Calspan  version 
of  the  MSCAM  nodel;  data  describing  the  number  of  bases  and  depots  in  the 
problea  and  the  number  of  aissiles  at  each  base;  and  data  defining  the  con- 
straint for  total  systea  LCC  (unscheduled  Maintenance). 

The  ARIES  aodel  processes  the  individual  SRU  data  independently  to 
derive  the  optiaua  stockage  levels,  Maintenance  strategy,  and  expected  back- 
order sua  across  all  bases  as  a function  of  the  LCC  (unscheduled  Maintenance) 
for  the  iten.  The  resultant  data  froa  the  individual  SRU  calculations  are 
stored  on  a penanent  disc  file  for  later  processing  in  the  Model.  The 
peraanent  file  eliainates  redundant  SRU  computations  and  provides  a data  bank 
for  future  analyses. 

Using  a resource  allocation  algoritha,  ARIES  confutes  the  optimal 
logistics  configuration  for  each  SRU  that  will  provide  the  greatest  amount 
of  systea  availability  relative  to  an  input  value  for  total  LCC  (unscheduled 
Maintenance).  The  solution  at  the  cost  constraint  defines  the  optiaua  Main- 
tenance strategy,  stockage  levels,  and  expected  backorder  value  for  each  SRU. 
By  supplying  the  aodel  with  the  proper  execution  option,  useful  suaaary  data 
nay  be  computed  for  each  base  and  also  at  the  depot (s). 

Since  the  LCC  for  logistics  is  only  one  part  of  the  total  LCC  fOT  a 
weapon  systea,  the  results  froa  an  analysis  of  istscheduled  maintenance  should 
be  aerged  with  the  costs  for  other  phases  of  the  systea  to  influence  the  total 
LCC  to  be  optimized.  To  illustrate,  there  is  a point  beyond  which  increasing 
the  systea  acquisition  cost  of  an  SRU  to  improve  reliability  and  Maintainabil- 
ity characteristics  will  not  be  offset  by  a corresponding  decrease  in  the 
unscheduled  Maintenance  cost. 


*There  are  four  alternative  maintenance  strategies:  base-depot  repair, 

depot-only  repair,  base-only  repair,  and  discard  ipon  failure. 


While  ARIES  was  developed  for  use  in  analyzing  the  logistics  of  a 
weapon  systen,  the  model  is  also  applicable  to  other  non-defense  logistics 
systems  that  have  a similar  base-depot  maintenance  concept  involving  high-cost, 
low-demand  recoverable  parts  (e.g.,  a computer  system  network  or  a telecommuni- 
cations system). 

USE  OF  THE  EXPECTED  NORS  VS  LCC  (UNSCHEDULED  MAINTENANCE)  CURVE 

An  important  curve  may  be  constructed  that  shows  the  relationship 
between  optimum  system  availability  (logistics)  and  LCC  (unscheduled  maintenance). 
The  curve  is  derived  by  solving  the  logistics  allocation  problem  with  ARIES  for 
a range  of  cost  constraints  that  cover  the  region  of  feasible  solutions.  Figure 
0 illustrates  the  principal  curves  that  were  constructed  from  the  analysis  of 
the  example  problem,  which  is  discussed  in  more  detail  in  Section  IV  and 
Supplement  A.  The  curves  of  Figure  0 show  how  the  expected  number  of  NORS 
missiles,  for  both  the  complete-cannibalization  and  no-cannibalization  cases, 
decreases  as  the  expenditure  for  LCC  (unscheduled  maintenance)  increases.  The 
dashed  line  in  the  figure  represents  the  expected  backorder  curve  until  it 
essentially  coincides  with  the  expected  NORS  curve  for  the  no-cannibalization 
case.  In  the  Air  Force  logistics  planning  activity,  reliance  is  placed 
primarily  upon  the  expected  NORS  for  the  case  where  no  cannibalization  is 
assumed. 

In  using  Figure  0,  the  logistician  should  select  a design  point  on 
the  curve  that  at  least  meets  the  weapon  system's  specification  for  logistics 
availability  and  that  also  takes  into  consideration  the  marginal  return  from 
additional  logistics  investment  (i.e.,  the  marginal  expenditure  for  LCC  to 
reduce  the  expected  NORS  by  one  does  not  exceed  the  life  cycle  cost  to  procure 
an  additional  missile).  If  one  were  to  consider  "some"  cannibalization,  the 
curve  for  it  would  lie  inside  the  two  curves  for  the  cases  assuming  no 
cannibalization  and  complete  cannibalization. 
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WEAPON  SYSTEM  LCC  (UNSCHEDULED  MAINTENANCE),  t MILLION 
Fig.  0.  ARIES  Anetysn  of  System  Availability  vt  LCC  for  Unecheduled  Maintananca 


USE  OF  THE  ITEM  EXPECTED  BACKORDER  VALUES 


The  planner  has  in  effect  selected  a LCC  (unscheduled  awintenance) 
constraint  point  after  he  has  perforaed  the  analyses  discussed  above.  Corre- 
sponding to  the  cost  constraint  point,  there  is  a solution  that  defines  the 
expected  backorder  value  for  each  SRU.  The  SRU  expected  backorder  value  is  a 
necessary  input  to  the  MSCAM  nodel  to  perfona  refined  reliability/naintainability 
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trade-off  analyses  on  each  individual  item.  The  expected  backorder  values  for 
the  sixteen  SRUs  in  the  example  problem  ranged  from  0.042  to  1.476  at  a 
cost  constraint  of  $106  M. 

MANPOWER  REQUIREMENTS 

The  ARIES  model  has  an  option  to  compute  the  unscheduled  maintenance 
manpower  requirements  at  each  base  within  each  skill  category  for  each  SRU. 
Also,  the  depot  manpower  requirements  are  computed  foT  each  SRU.  The  total 
base  manpower  requirement  and  the  total  depot  manpower  requirement  are  also 
summarized. 

COMMON  AGE  (AEROSPACE  GROUND  EQUIPMENT) 

A limitation  with  both  MSCAM  and  ARIES  is  that  items  are  treated 
independently;  hence,  the  models  do  not  have  explicit  logic  to  take  into 
account  the  fact  that  several  items  may  use  common  AGE  for  unscheduled  main* 
tenance.  Three  approaches  are  identified  below  to  assist  in  remedying  this 
situation. 

Grouping  items  - It  may  be  possible  to  group  several  items  that 
use  the  same  AGE  and  have  the  same  maintenance  strategy  and 
other  important  characteristics  and,  thus,  to  treat  these 
items  in  the  model  as  if  they  were  one  item. 

Prorating  AGE  expenses  - It  may  be  possible  to  estimate  the  pro- 
portionate share  of  common  AGE  usage  by  several  items.  The 
input  data  to  the  model  would  reflect  the  pro  rata  AGE  expenses. 

Attributing  entire  AGE  expense  to  one  item  - It  may  be  appro- 
priate to  attribute  the  entire  AGE  expense  to  only  one  item 
from  among  several  that  potentially  could  use  the  same 
common  AGE.  This  situation  may  be  justifiable  if  one  '.tern 
has  been  planned  for  AGE-assisted  repair  in  any  eventuality. 

It  may  be  necessary  to  iterate  through  the  solution  several  times  to  arrive  at 
the  proper  procedure  for  accomaodating  the  common  AGE  expenses. 
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TYPES  OF  TRADE-OFFS  THAT  CAN  BE  MADE 


Obviously,  the  types  of  trade-offs  that  can  be  made  with  any  model 
are  variations  of  the  many  combinations  of  input  variables.  For  the  ARIES 
model,  it  is  convenient  to  summarize  the  principal  parameters  for  trade-off 
analyses  according  to  the  following  categorization: 

Reliability 

e Failure  rates 

• Operating  hours  per  base  per  month 
Maintainabi lity 

e Base  and  depot  AGE 

• Base  and  depot  repair  capability: 

(1)  Man-hours  per  repair 

(2)  Repair  parts  cost  per  repair 

(3)  Repair  turnaround  time 

Availability 

e System  deployment  concept 
e Expected  item  backorders  across  all  bases 

e Expected  NORS  vehicles 

Item  Charaoterietioe 

• End-item  cost 

e NRTS  (Not  Repairable  This  Station)  and  condemnation  rates. 


Section  nr 

TECHNICAL  DISCUSSION 


ARIES  is  a mathematical  model  for  the  multi-echelon  base- 
depot  supply  system  for  recoverable  items.  Item  demand  is 
Poisson.  There  aTe  four  alternative  maintenance  strategies 
for  an  item  and  the  strategy  choice  is  a variable  in  the  prob- 
lem. Rather  than  considering  only  an  item's  unit  cost,  the 
model  computes  the  item's  10-year  life  cycle  cost  for  unsched- 
uled maintenance.  The  sum  of  item  expected  backorders  across 
all  bases  is  minimized  subject  to  a life  cycle  cost  constraint 
for  unscheduled  maintenance.  The  item  expected  backorder  versus 
cost  function  can  have  large  regions  of  non-convexity.  The 
resource  expenditure  algorithm  for  the  model  incorporates 
allocations  from  these  non-convex  regions  to  provide  solu- 
tions near  a constraint.  The  maximum  region  of  uncertainty 
for  the  minimized  expected  backorder  sum  is  defined  and  this 
provides  a basis  for  judging  the  quality  of  the  result. 

ARIES  evaluates  the  minimized  expected  backorder  solution  to 
determine  the  expected  number  of  NORS  vehicles  for  both  the 
case  where  coiplete  cannibalization  of  parts  is  assumed  and 
the  case  where  there  is  no  cannibalization.  The  manpower 
and  item  stockage  level  requirements  for  the  individual  bases 
and  depot  are  determined  by  the  ARIES  solution. 
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THIS  PAPER  discusses  • new  logistics  Model,  ARIES,  that  combines  and 
extends  the  features  of  two  previously-developed  Models  to  achieve  a More 
general  analytic  treatment  for  the  multi-iten,  multi-echelon  supply  system 
for  recoverable  parts.  The  M2TRIC  model  ^ analyzed  aggregate  recoverable 
item  problems  while  considering  only  an  item's  unit  cost  and  parameters  that 
defined  a single  maintenance  strategy  for  the  item  for  all  ranges  of  resource 
allocations.  The  MS  CAM  (Manpower,  System  support  Cost  Analysis  Model)  model 
has  as  its  objective  the  analysis  of  an  individual  recoverable  item  for  each 
of  four  alternative  maintenance  strategies  to  determine  the  10-year  life  cycle 
cost  (unscheduled  maintenance)  and  the  base  and  depot  stockage  levels  corres- 
ponding to  a given  input  value  for  the  item  sum  of  expected  backorders  across 
all  bases. 

The  military  services  have  been  under  increasing  pressure  in  recent 
years  to  estimate  and  reduce  the  cost  of  developing,  deploying,  and  opera- 
ting new  sophisticated  weapon  systems.  The  "cost  growth"  phenomenon  of  past 
weapon  systems  has  brought  about  enough  attention  such  that,  once  the  require- 
ments have  been  established  for  a new  weapon  system,  the  following  order  of 
priorities  tends  to  prevail:  (1)  cost,  (2)  performance,  and  (3)  schedule. 

The  logistic  support  costs  for  a weapon  system  account  for  a signifi- 
cant portion  of  the  total  weapon  system  cost  over  the  life  cycle  of  the  deploy- 
ed system.  In  1964,  the  Department  of  Defense  issued  Directive  4100.35  des- 
cribing the  Integrated  Logistics  Support  (ILS)  concept.  Since  then  the  Services 
have  been  required  to  consider,  estimate  and  evaluate  the  life  cycle  costs 
(LCC)  Implied  by  the  design  alternatives  encoisitered  throughout  the  acquisi- 
tion process.  By  providing  improved  methods  for  making  trade-off  evaluations 
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early  in  the  weapon  systems  design  phase,  attention  can  be  focused  upon 
reducing  life  cycle  costs  before  significant  funds  are  connitted  to  a 
design  or  concept. 
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Within  the  SCAD  (Subsonic  Cruise  Armed  Decoy  for  the  B-52  aircraft) 
program,  the  MSCAM  model  was  modified  for  use  to  perform  refined  reliability/ 
maintainability  parameter  trade-offs  for  individual  recoverable  items  that 
are  called  SRUs  (Shop  Replaceable  Units) . Since  a key  input  variable  to  the 
MSCAM  model  is  the  expected  backorder  value  for  each  item,  the  requirement 
existed  to  develop  an  aggregate  model  to  derive  the  target  value  for  the  item 
sum  of  expected  backorders  across  all  bases  when  one  considers  LCC  and  four 
alternative  maintenance  strategies.  Another  feature  the  aggregate  model  re- 
quired is  that  it  relate  to  an  operational  measure  of  effectiveness  for 
system  availability. 

Since  the  LCC  for  logistics  is  only  one  part  of  the  total  LCC  for  a 
weapon  system,  the  results  iron,  an  analysis  of  unscheduled  maintenance  should 
be  merged  with  the  costs  for  other  phases  of  the  system  to  influence  the  total 
LCC  to  be  optimized.  To  illustrate,  there  is  a point  beyond  which  increasing 
the  system  acquisition  cost  of  an  SRU  to  improve  reliability  and  maintainability 
characteristics  will  not  be  offset  by  a corresponding  decrease  in  the  unscheduled 
maintenance  cost. 
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While  ARIES  was  developed  for  use  in  analyzing  the  logistics  of  a 
weapon  system,  the  model  would  be  applicable  to  other  non-defense  logistics 
systems  that  have  a similar  base-depot  maintenance  concept  involving  high-cost , 
low-demand  recoverable  parts,  e.g.,  a computer  system  network  or  a telecommuni- 
cations system. 
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MODEL  OBJECTIVES 


ARIES  is  a computerized  mathematical  model  to  assist  the  ILS  manager 
in  the  decision-making  process  for  optimizing  system  component  reliabili- 
ties, maintenance  strategies,  and  deployment  stockage  levels  while  taking 
into  account  the  corresponding  impact  upon  system  availability  and  LCC  for 
unscheduled  maintenance.  For  brevity,  the  term  LCC  shall  refer  to  life  cycle 
cost  for  unscheduled  maintenance  for  the  remainder  of  the  paper.  Subject  to 
a given  cost  constraint,  ARIES  minimizes  the  expected  sum  of  item  backorders 
across  all  bases.  The  operational  measure  of  effectiveness  to  which  availa- 
bility is  best  related  to  by  operational  personnel  is  the  expected  number  of 
NORS  (Not  Operationally  Ready  - Supply)  vehicles.  The  solution  for  the  expect- 
ed backorder  result  can  be  evaluated  to  determine  the  corresponding  expected 
number  of  NORS  vehicles.  By  the  systematic  evaluation  of  a range  of  cost 
constraints,  a useful  curve  may  be  derived  to  show  the  relationship  between 
availability  and  LCC. 

A fallout  of  the  LCC  analysis  is  the  determination  of  the  number  of 
men  required  in  each  skill  category  at  the  bases  and  depot  to  perform  un- 
scheduled maintenance.  ARIES  is  structured  as  a computer  model  that  realizes 
the  benefits  from  creating  a permanent  data  file  for  individual  item  expect- 
ed backorder-LCC  functions  and  this  facilitates  future  analyses  by  eliminating 
the  need  to  recalculate  all  the  item  data  when  the  same  basic  system  is  under 
study  md  only  a few  parameters  are  being  changed. 

MODEL  ASSUMPTIONS 

The  mathematical  assumptions  describing  the  structure  of  the  model  are 
given  below.  Since  many  of  the  assumptions  are  the  same  as  those  previously 
described  for  METRIC!8!,  only  a brief  mention  of  them  will  be  made  here  to  de- 
fine the  problem. 
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Item  Demand.  Ie  Poisson 


The  probability  distribution  characterizing  item  demand  arrivals  is 
assumed  to  be  Poisson.  Considering  the  typical  accuracy  of  item  failure 
data  that  is  available,  particularly  during  the  design  phase  of  a weapon 
system,  the  simple  Poisson  distribution  represents  a reasonable  choice  for 
the  model  (also,  see  reference  9,  p.  632).  As  a consequence  of  this  assump- 
tion, demand  is  a stochastic  process  with  independent  increments. 

Multi-Eohelon  Repair 

Item  demand  may  be  satisfied  at  either  a base  or  depot.  Each  item 
has  one  depot.  The  depot  need  not  be  the  same  for  all  items.  The  mainten- 
ance strategy  that  is  selected  for  an  item  will  determine  where  repair  is 
permitted. 

Four  Alternative  Maintenance  Strategies 

Each  item  will  have  one  of  the  following  maintenance  strategies  corres- 
ponding to  a given  resource  expenditure  level:  repair  at  both  base  and  depot, 

(BD) ; repair  at  the  depot  only,  (D) ; repair  at  the  base  only,  (B) ; and 
no  repair-discard  upon  failure,  (NR) . 

When  a demand  occurs  at  the  base,  it  will  be  repaired  there  with  proba- 
bility 1-r  or  sent  to  the  depot  for  repair/replacement  with  probability  r. 

In  the  vernacular  of  logistics,  r is  referred  to  as  the  NRTS  rate  (Not 
Reparable  This  Station).  For  the  D and  NR  strategies,  the  value  of  r is 
one.  When  applying  the  multi -echelon  theory  to  compute  an  item's  expected 
backorder  sum,  it  is  necessary  to  estimate  the  average  depot  resupply  time. 

The  depot  resupply  time  is  the  sum  of  the  average  delay  for  depot  repair  and 
the  round  trip  base -depot -base  order  and  shipping  time.  For  the  B and  NR 
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strategies,  the  "depot  repair  tine"  is  equal  to  the  average  factory  re- 
order tine  for  the  item. 

(8-1,8)  Policy  at  Base 

The  one-for-one  replacement  inventory  policy,  (s-l,s),  is  assumed  at 
the  bases. 

Condemnations 

For  the  B and  NR  maintenance  strategies  the  condemnation  rates  are  r 
and  1.0,  respectively,  and,  hence,  base  condemnations  and  subsequent  depot 
replacements  are  treated  by  the  (s-l,s)  policy.  The  condemnation  rate  at 
the  depot  for  the  BD  and  D strategies  is  an  input  parameter  that  is  usually 
under  5 percent.  The  model  considers  the  effect  of  condemnations  for  the 
BD  and  D strategies  only  to  the  extent  that  it  affects  recurring  stockage 
costs  and  ignores  the  small  effect  it  would  have  on  the  depot  Tesupply  time. 

No  Lateral  Resupply 

Simplifying  the  problem  is  the  assumption  that  there  is  no  lateral 
resupply  among  the  bases.  During  the  planning  phase,  it  is  not  Air  Force 
policy  to  consider  the  lateral  exchange  of  stock  as  a source  of  inventory. 

Item  Demand  Ie  Stationary  During  the  Life  Cycle  Period 

While  this  assumption  may  be  more  difficult  to  accommodate  for  air- 
craft systems  that  have  a widely  varying  use  rate,  it  is  quite  appro- 
priate for  many  missile  systems,  since  the  demand  rate  during  the  life 
cycle  period  is  relatively  constant. 

Reeupply  Time  Ie  Independent  of  Demand 

The  time  to  resupply  an  item  at  the  base  or  depot  is  independent 
of  the  demand  occurring  (infinite  channel  queueing  model) . 


17 


r 


r 


s 


I 


Baae-To-Depot  Item  Demand  Can  Be  Pooled 

The  net  effect  of  the  demand  from  all  the  bases  to  the  depot  can  be 
pooled  to  form  a resultant  new  Poisson  distribution.  The  depot  demand  rate 
is  the  sum  of  the  individual  base-to-depot  demand  rates. 

Equal  Baee  and  Item  Eaaentialitiee 

ARIES  is  presently  structured  to  treat  each  base  and  each  item  with 
the  same  importance.  Any  item  can  cause  a NORS. 

Minimum  of  One  Unit  of  Stock  Per  Item 

The  model  assumes  the  system  has  a minimum  of  at  least  one  unit  of 
stock  per  item  in  inventory.  This  assumption  facilitates  confuting  the 
minimum  LCC  for  the  system  and  the  corresponding  expected  NORS.  Logistics 
planners  indicate  that  operational  managers  would  require  there  be  at  least 
one  unit  per  item  in  inventory. 

Minimizing  the  Sum  of  Baee  Expected  Backordere 

A demand  that  can  not  be  satisfied  from  base  supply  is  a backorder 
or  due-out.  The  backorder  measurement  takes  into  account  the  duration  of 
a demand  as  well  as  the  fact  that  it  occurred.  The  model  minimizes  the 
sum  of  item  expected  backorders  across  all  bases.  As  will  be  discussed 
later,  the  process  of  minimizing  the  sum  of  expected  backorders  provides 
an  excellent  approximation  to  the  process  of  minimizing  expected  NORS 
for  the  case  when  there  is  no  cannibalization. 

INPUT  DATA 

The  input  data  to  the  model  can  be  categorized  into  data  that  is 
common  to  all  items,  individual  item  data,  and  peculiar  problem  data.  The 
following  inputs  are  included  in  ARIES.  Many  of  the  parameters  are 
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necessary  to  compute  the  non-recurring  and  recurring  costs  that  go  into 
the  LCC  computations. 

Data  Caiman  To  All  Item 

Total  operating  hours  for  all  vehicles  at  eaoh  base  per  month. 
Base-depot  distances  and  order  and  shipping  time. 

Base  and  depot  manhour  ooet  rates. 

Packaging  and  shipping  rate  oosts. 

Ten-year  uniform  series  present  worth  factor.  This  factor  is  neces 
sary  for  combining  the  annual  recurring  cost  with  the  non-recurring  cost 
to  determine  the  present  value  for  the  10-year  LCC. 

Item  Data 


Failure  rate.  The  equipment  "on"  and  "off"  time  failure  rates  are 
reduced  to  one  number  that  when  multiplied  by  the  total  vehicle  operating 
hours  and  the  number  of  item  units  per  vehicle  yields  the  expected  item 
demand  per  base  per  month. 

Base  and  depot  AGE  (Aerospace  Ground  Equipment)  cost  per  unit. 

AGE  installation  and  set-up  cost. 

Annual  operating  and  maintenance  cost  per  AGE. 

Average  number  of  hours  that  depot  AGE  is  used  per  repair. 

Item  unit  cost,  shipping  weight,  EFTS  rate , and  condemnation  rate. 

Base  and  depot  average  repair  time  manhours. 

Base  and  depot  average  repair  turnaround  time. 

Base  end  depot  repair  parts  cost  per  repair.  This  is  the  cost  of 
consumable  parts  per  repair. 

Average  base  manhours  to  detect,  to  remove,  and  to  replace  a failed  item. 
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Proportion  of  repairs  performed  by  each  of  the  baee  labor  ekill  categories. 

Nunber  of  pages  of  technical  data. 

Procurement  lead  time. 

Nunber  of  units  of  the  item  used  per  vehicle. 

Problem  Data 

Lagrange  multipliers. 

LCC  and  sum  of  expected  backorder  constraints. 

Nunber  of  vehicles  at  each  base. 

A discussion  of  the  calculations  that  go  into  co^niting  the  LCC  is 
found  in  reference  S.  Briefly,  the  non-recurring  and  recurring  costs  are  a 
summation  of  the  costs  for  stockage,  AGE,  transportation,  and  repair  (labor 
and  materials) . The  ARIES  model  has  provision  to  accommodate  a fraction  of 
the  BD  and  B maintenance  strategy  base  repair  and  AGE  costs  for  the  D and 
NR  strategies.  This  feature  reflects  the  cost  at  the  base  to  detect,  to 
remove,  and  to  replace  failured  items  that  are  either  sent  to  the  depot 
or  discarded. 

ANALYTIC  STRUCTURE  OF  ARIES 

The  model  has  three  principal  logical  steps.  The  first  step  performs 
a sub -optimization  to  derive  the  individual  item  expected  backorder-LCC 
functions.  The  second  step  employs  Lagrange  multiplier  and  marginal  allo- 
cation techniques  to  reach  a solution  at  a problem  constraint.  And  last, 
the  third  step  evaluates  the  expected  backorder  solution  at  a constraint  in 
terms  of  its  corresponding  expected  NORS  values. 

Mathematical  Problem  Definition 

ARIES  solves  either  of  two  single-constraint  problems  that  require 


It  has  been  shown1  J that  problems  (1)  and  (2)  can  be  refoimulated 
into  the  fora  of  equation  (3)  when  the  item  expected  backorder-LCC  func- 
tions are  convex.  In  the  case  of  ARIES,  we  are  using  the  item  LCC  instead 
of  the  item  unit  cost,  but  this  in  no  way  affects  the  validity  of  equation 
(3).  Thus,  find  (s^}  «»d  (■j)  which  minimizes 

F - Y.  B..  ♦ LCC  , (3) 

i» j 13  i 

where  X > 0 is  an  unknown  Lagrange  multiplier.  Since  F can  be  treated 
as  a separable  cell  integer  programming  problem,  one  can  take  advantage 
of  restricting  attention  to  one  item  at  a time:  minimize 

F.  - Y Bjj  ♦ XLCCj,  (i  » 1,2,...  ,1) . (4) 

Summing  the  Fi#  F ■ ^Fj.  The  multi-echelon  theory  that  derives  the  basic 
expression  for  is  presented  in  [e]  • A theorem,  useful  for  conqputational 
purposes,  on  the  recursive  expression  for  an  item  expected  backorder 
sum  is  given  in  the  appendix. 

Item  Expected  Baokorder-LCC  Function 

To  facilitate  solving  equation  (4),  it  is  first  necessary  to  calculate 
for  each  item  and  for  each  of  the  four  maintenance  strategies  the  set  of 
allocations  which  are  undominated  - in  the  sense  that  each  point  ( 6.  ., 

i 3 

LCC)  of  that  set  in  the  expected  backorder— LCC  plane  is  either  below  or  to 
the  left  of  every  other  allocation  point.  This  set  is  called  the  item  expected 
backorder — LCC  function  (pure  maintenance  strategy).  Figure  1 illustrates 
how  the  four  pure  maintenance  strategy  expected  backorder  versus  LCC 
discrete  functions  could  look  for  an  item.  The  functions  are  monotonically 
decreasing  and,  in  this  example,  they  alternately  dominate  each  other  at 
different  regions  over  the  range  of  LCC. 


■ 4 


I 

I 

I 
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MAINTENANCE 

STRATEGIES 


Fig.  1. 


LCC.  S 


Sum  of  Expected  Back  or  dan  vs  LCC  for  Each  of  Four  Altamativa  Maintenance 
> for  Item  (i) 


I 

I 

I 

I 
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The  aodel  collapses  the  four  functions  into  one  suboptinal  expected 
backorder  sua-LCC  function  for  ite»  i;  call  it  {Cjn)  where  Cjn  represents 

ZeU 

Bj  j , LCCj)  at  the  n point  in  the 

undoainated  set  in  the  direction  of  increasing  LCC.  The  circled  points 
in  Figure  1 illustrate  the  undoainated  set  in  the  exaaple.  Within  the 
coaputer  program  for  the  aodel,  a large  savings  in  coagmtation  is  realized 
by  testing  whether  a previously-calculated  strategy  doainates  the  one 
under  current  processing.  ARIES  computes  the  strategies  in  the  order  of 
BO,  D,  B and  NR  and  this  is  generally  the  order  of  decreasing  selection 
frequency . 
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Convexifioation  of  the  Item  Expected  Backorder-LCC  Function 

Unfortunately,  it  turns  out  that  the  item  expected  backorder-LCC 
function  can  have  large  regions  of  non-convexity.  Three  potential  sources 
of  non-convexity  arise  from  the  allocation  of  an  additional  unit  of  total 
stock  for  the  item  causing  (1)  a shift  in  the  amount  of  stock  at  the  depot, 
(2)  a change  in  maintenance  strategy,  or  (3)  a disportionate  reduction  in 
expected  backorders  because  of  unequal  base  demand  rates. 


The  item  expected  backorder-LCC  function  is  modified  to  form 

a new  function  (cjp)  that  includes  only  those  of  the  originally  defined 
points  that  are  on  the  lower  left  boundaxy  of  the  convex  hull  of  {CjnL 
The  new  consecutive  sequence  of  points  of  increasing  LCC  are  indexed  by 
the  subscript  p. 


Iteration  To  a Problem  Constraint 

In  the  second  principal  logical  step  there  are  three  stages  in  the  process 
of  minimization  under  a problem  constraint.  For  a given  Lagrange  multiplier  X 
and  the  convex  functions  {cjp),  (i  ■ 1,2..., I),  equations  (3)  and  (4)  may  be 
solved,  i.e.,  minimize  the  and,  thus*  F.  This  solution  is  representable 
by  the  I-  tuple  of  integers  (pj,  p2,  ....  Pj)  identifying  the  selected  points 
from  each  item  expected  backorder-LCC  function,  i.e.,  ( ( 51  Bjj)p  » (LCCi)p.)# 
If  equations  (3)  and  (4)  are  solved  for  a given  set  of  Xfc,  (k  ■ 1,2, . . . ,K) , then 
a test  can  be  made  to  determine  between  which  two  X^  a problem  constraint  is 
met;  call  the  two  bounding  Lagrange  multipliers  Xml  and  Xn< 
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The  model  employs  incremental  marginal  allocation  with  the  functions 
(Cjp)  starting  from  the  point  in  resource  expenditure  defined  by  the  solu- 
tion at  The  marginal  allocation  procedure  iteratively  selects  and 

allocates  to  the  item  i for  which  the  quantity  ^ , as  given  by  equation 
(5), is  a maximum  and  stops  when  a constraint  would  be  first  violated. 

♦i  - - ( Z f^ql/iaccpq  - CLCCp^j],  where  q-Pj.  (5) 

Whenever  an  allocation  is  made,  q » p.  is  increased  by  one.  The  solu- 
tion generated  by  the  marginal  allocation  with  ^ is  an  optimal  allo- 
cation of  resources  for  its  corresponding  expenditure  (or  expected  back- 
order sum);  however,  the  solution  is  not,  in  general,  optimal  for  the 
problem  constraint. 

Starting  with  the  solution  determined  by  marginal  allocation  with  ^ 

from  equation  (5) , a heuristic  marginal  allocation  algorithm  is  next 

applied,  but  now  using  all  points  from  the  original  functions  {<;.  } . Hie 

in 

procedure  iteratively  selects  and  allocates  to  the  item  i for  which  the 
quantity  e^,  as  given  by  equation  (6),  is  a maximum  until  there  is  no 
remaining  item  allocation  left  that  will  not  violate  a constraint. 


(6) 


“ [(  Z ' ( Z 8. j)t]/[(LCCi)t  - (LCC1)t_1]»  where  t-n^ 

Whenever  an  allocation  is  made,  t«n^  is  increased  by  one.  The  I-tuple  of 
integers  (n^.nj,. . . ,nj)  identifies  the  solution  at  any  point  in  the  final 
allocation  process. 

The  final  marginal  allocation  procedure  does  not  guarantee  that  the  exact 
optimal  solution  will  be  found  since  allocations  may  be  made  from  non- convex 
regions  of  the  (c^)  functions.  The  maximum  region  of  uncertainty  for  the  mini- 
mized expected  backorder  sum  can  be  defined,  and  this  provides  a basis  for 
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*1  *u 


evaluating  the  quality  of  a resulting  solution.  Using  the  functions  and 

a sufficiently  large  set  of  the  entire  array  of  Lagrange  multiplier  solu- 
tions may  be  derived.  The  circled  points  in  Figure  2 represent  how  the  solu- 
tions would  appear.  These  solutions  have  the  following  properties: 


(i)  Connecting  adjacent  solutions  by  straight  line  segments  yields 
a convex  polygonal  line;  and 

(ii)  No  solution  corresponding  to  either  the  convexified  functions 

or  the  original  item  functions  Uin>  can  occur  below  the  connecting 
convex  polygonal  line. 


I tec, 

Maximum  Aran*  of  Uncertainty  for  Optimal  Solutions  Using  the  Original  Item 
Expected  Oechorder  - LCC  Functions 


Property  (i)  follows  from  the  fact  that  when  the  Lagrange  multiplier 
fonulation  is  applied  incrementally  (iterative  marginal  allocation)  to  a 


set  of  convex  functions  the  resultant  solutions  constitute  a new  function 
that  is  also  convex,  i.e.,  the  function  and  slope  will  be  monotonically 
decreasing  and  increasing,  respectively.  Property  (ii)  is  true  for  the 
(c!p)  functions  by  the  very  nature  of  Lagrangian  analysis  providing  an  un- 
dominated sequence  of  solutions.  Property  (ii)  is  true  for  the  func 

tions  since  allocations  from  non-convex  regions  must  lie  above  and  to  the 
right  of  the  tangent  lower  bounding  hyperplanes  (convex  polygonal  line  in 
ARIES)  to  the  envelope  for  the  entire  set  of  all  possible  solutions  in  the 


By  connecting  the  Lagrangian  and  the  final  heuristic  marginal  alloca- 
tion solutions  with  horizontal  and  vertical  lines,  as  shown  in  Figure  2,  the 
area  enclosed  between  these  lines  and  the  convex  polygonal  line  defines  the 
region  of  maximum  uncertainty  for  optimal  solutions  when  considering  the 
original  item  expected  backorder-LCC  functions  • The  maximum  region 

for  uncertainty  for  the  solution  corresponding  to  a cost  constraint  is 
illustrated  in  Figure  2 by  the  cross-hatched  area.  The  difference  between 
the  solution  that  is  obtained  and  the  theoretical  minimum  possible  expected 
backorder  sum  is  denoted  by  AB  in  the  Figure. 

Expected  NORS  Evaluation 

The  third  principal  logical  step  in  the  model  takes  the  item  stockage 

* 

levels  corresponding  to  a problem  constraint  minimized  expected  backorder 
sum  solution  and  computes  the  expected  number  of  NORS  vehicles  for  both  the 
case  when  complete  cannibalization  is  assumed  and  the  case  where  there  is 
no  cannibalization.  The  sum  of  expected  NORS  for  the  case  with  complete 


cannibalization  is  a function  of  the  item  in  shortest  supply  at  each  of 
the  bases. 


For  each  base  j,  define  ip^  by  the  following  expression: 


*j  * maxfg^/u^, 
i « 1 < i 


(7) 


where  is  the  number  of  units  of  item  i used  per  vehicle.  The  item 
in  shortest  supply  at  base  j is  defined  as  that  i which  maximizes  the 
expression  in  equation  (7).  Summing  ip.  over  all  the  bases  yields  the  ex- 
pected NORS  for  the  cannibalization  case 

E(NORS) 


?V 


(8) 


Contrary  to  Sherbrooke's  tentative  conclusion^ , the  process  of  minimizing 

the  sum  of  expected  backorders  provides  only  a fair  approximation  to  the 

process  of  minimizing  EfNORS)  . Attacking  the  minimization  of  E(NORS)  more 

^ c 

directly  with,  a final  allocation  algorithm  that  identifies  the  next  expendi- 
ture by  selecting  the  item  that  has  a maximum  for  the  quantity  &. ./u.  , 

j 13  1 

tends  to  generate  cost  constraint  solutions  that  have  a significantly  lower 


one 


value  for  EfNORS)  . However,  the  corresponding  expected  backorder  sum  is 

significantly  larger  than  that  which  is  obtained  by  the  marginal  allocation 

algorithm,  i.e.,  equations  (5)  and  (6).  Miller^  has  shown  that  the  EfNORS)  - 

c 

cost  function  is  convex  when  restricting  the  solution  to  fc^l  and  that  any 
one  of  several  non-linear  programming  techniques  could  be  applied  to  find  the 
optimum.  The  example  discussed  in  a later  section  illustrates  the  extent  of 
the  EfNORS) c differences  between  the  two  ARIES  algorithms  for  approaching  a 
problem  constraint.  Since  in  the  Air  Force  logistics  planning  activity  reli- 
ance is  placed  primarily  upon  the  expected  NORS  for  the  case  where  no  canni- 
balization is  assumed,  EfNORS)  , it  was  considered  sufficient  and  practical 

IlC 


^«s®*ws?5 
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to  use  the  algorithm  that  employes  £ B../u.  in  the  model  as  an  approxi- 

i 3 

mat ion  for  the  minimum  value  of  E(NORS)c. 

Underlying  the  Air  Force  reliance  on  E(NORS)  is  the  concept  that  can- 

nc 

nibalization  should  only  be  used  as  a fall-back  technique  in  the  event  that 
actual  item  demand  and/or  resupply  times  are  higher  than  previously  anti- 
ciapted.  The  value  for  E(NORS)nc  is  derived  as  follows: 

Let  Pj  » Pr  (the  first  unit  selected  for  item  i on  a random  vehicle 
at  base  j is  backordered) , 

P2  * Pr  (the  second  unit  selected  for  item  i on  a random  vehicle 
at  base  j is  backordered | first  unit  is  not  backordered). 


P * Pr  (the  u.  unit  for  item  i on  a random  vehicle  at  base  j 
u^  v 1 

is  backordered | the  previous  u^-1  units  are  not  backordered). 


Vj  » the  nunber  of  vehicles  at  base  j,  then 

P1  * 6i j/vjui * P2  " 6ij/(vjui“1)  * •••*  pUi  " Bi/V/lVYl)*!].  (9) 
where  is  the  expected  backorders  for  item  i at  base  j when  the  demand  rate 
is  decreased  by  the  fact'-  (VjU^-k+l)/v^u^ . Assuming  independence  of  item 
unit  demand,  the  probabilit.  it  a vehicle  is  not  NORS  due  to  item  i at  base  j is 


* ■ K>v- 


Since  it  would  require  extensive  computations  to  calculate  ?2»  P3»  • • • » pUi * 
and  the  decrease  in  the  numerator  will  be  approximately  equal  to  the  decrease 
in  the  denominator  in  equation  (9) , the  simplifying  approximation  is  made 
that  ■ p2  ■ ...  - pUi-  Assuming  item  demand  independence,  the  probability 


j is  Tf  A.  . 
i « 

7T  A.,  a TT  (1-6^/v  u.)ui 
i J i J 


a vehicle  is  not  NORS  at  base  j is  II  A.,  where, 
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(11) 


Hence,  the  expected  total  number  of  NORS  vehicles  across  all  bases  is 
given  by  equation  (12) . 


E(NORS)nc  s Z v.(l-  IT  A..) 


(12) 


The  E(NORS)nc  function  is  non-separable;  however,  the  process  of  mini- 
mizing the  sum  of  expected  backorders  provides  an  excellent  approximation 
to  the  process  of  minimizing  E(NORS)nc>  This  fact  can  be  shown  by  examining 
the  expansion  of  equation  (12)  and  rearranging  terms.  Initially,  assume 
all  Uj  ■ 1,  then 


E(NORS)nc 


I 

i -i 

I i_ 

I 6 

LJ 

1 yi 

i<k 

I 

j 

J 

i 

ij8kj  + ^v2.  £k</ij0kje‘j 


ij 


(13) 


Since  all  3.  < v.,  E(NORS)  is  equal  to  the  sum  of  expected  backorders  minus 

ij  J nc 

the  sum  of  terms  involving  increasing  orders  of  0^^ /v j up  to  the  order 
1-1  in  the  last  term.  The  effect  of  having  some  u^l  is  to  increase  the 
number  of  such  terms.  For  the  situation  when  B^j<<v^, 

E(N0RS1  3 Z (14) 

nc  ij 
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If  the  item  demand  rates  remain  unchanged  but  the  number  of  vehicles  Vj 
is  varied,  an  increase  in  v^  will  further  reduce  the  difference  between  the 
expected  backorder  sum  and  E(NORS)nc. 

EXAMPLE 

Consider  a problem  with  16  items,  17  bases,  one  depot,  all  u^  » 1, 
and  an  unequal  number  of  vehicles  at  the  bases.  The  other  input  values  to 
the  model  for  the  example  are  too  numerous  to  mention  here  but  let  it  be 
remarked  that  a representative  range  of  item  unit  costs,  demand  rates,  and 
AGE  costs  are  in  the  problem. 

The  results  for  the  solution  corresponding  to  eleven  cost  constraints 
are  summarized  in  Tables  I and  II.  The  first  cost  constraint  dictates  the 
solution  for  the  minimum  number  of  units  of  stock  for  the  system,  i.e.,  one 
per  item.  The  remaining  constraints  cover  a range  of  LCC  values  down  to 
the  point  where  there  is  very  little  improvement  left  that  can  be  made  in 
system  availability.  There  are  three  alternative  allocation  algorithms 
presented  in  Table  I and  they  are  defined  below: 

Algorithm  1 - Using  only  the  points  on  the  convexified  item  functions, 
, marginal  allocation,  i.e.,  equation  (5),  is  employed  until 
a constraint  would  be  first  violated. 

Algorithm  2 -Starting  from  the  algorithm  1 solution  and  using  the 
item  functions  {cjp>  • repeated  marginal  allocations  are  made 
until  there  remains  no  item  allocation  that  would  not  violate  a 
constraint . 


Algorithn  3 -Starting  from  the  algorithm  1 solution  and  using  the 
original  item  functions  (c^},  repeated  marginal  allocations  are 


TABLE  I 

COMPARISON  OF  ALLOCATION  ALGORITHMS  FOR  MINIMIZING 
THE  SUM  OF  EXPECTED  BACKORDERS 


COST 

CONSTRAINT 

t.M 

SUM  OF 

EXPECTED  BACKORDERS 

MAXIMUM  PERCENT  POOSIOLE 
REDUCTION  FOR  SOLUTION 

EINORS)^ 

i 

2 

3 

1 

2 

3 

1 

2 

_i 

§3.5 

162.0 

102.6 

102.6 

0.23 

0.23 

043 

1404 

§4 

112.1 

1114 

1114 

0.02 

0.12 

0.06 

106.0 

106.4 

96 

90.3 

09.1 

96.6 

6.63 

0.76 

0.06 

96.3 

81.0 

99 

07.0 

00.7 

064 

0.06 

0.12 

0.07 

644 

04.4 

99 

42.7 

42.6 

424 

0.27 

0.06 

0.06 

41.9 

414 

414 

100 

2M 

27.9 

27.7 

4.30 

0.76 

0.10 

20.4 

27.5 

274 

102 

17.7 

17.0 

17.6 

0.90 

0.14 

0.07 

17.6 

17.4 

104 

11.7 

11.0 

947 

140 

044 

114 

10.9 

100 

040 

6.02 

2.90 

0.27 

040 

6.70 

6.61 

100 

3.47 

3.47 

3.47 

0.07 

KB 

0.02 

3.47 

3.46 

110 

2.40 

240 

2.10 

1446 

KSfl 

046 

2.30 

2.19 

] 

TABLE  n 

COMPARISON  OF  ALLOCATION  ALGORITHMS  FOR  MINIMIZING  f 

THE  SUM  OF  EXPECTED  NORS  WHEN  CANNIBALIZING 


made  until  there  remains  no  item  allocation  that  would  not 
violate  a constraint. 


Algorithm  1 corresponds  to  the  best  Lagrange  multiplier  solution  that 

could  be  obtained.  Algorithm  3 is  the  basic  allocation  procedure  for  ARIES, 

i.e.,  equations  (S)  and  (6).  Algorithm  2 is  included  for  the  purpose  of 

examining  how  much  is  gained  by  using  the  non-convex  regions  of  the  item 

functions  {?.  } in  the  solution.  The  columns  in  Table  I for  "Maximum  Percent 
in 

Possible  Reduction  for  Solution"  are  related  to  A8  that  is  illustrated  in 
Figure  2.  For  this  example,  the  maximum  percent  possible  reduction  in  expected 
backorder  sum  for  algorithms  1,  2,  and  3 were  less  than  15,  5,  and  0.25  per- 
cent, respectively.  For  a different  problem  (that  is  not  presented  in  this 
paper)  that  had  many  items  with  u^  > 1 , the  corresponding  results  were  25, 

5,  and  1 percent  for  the  three  algorithms.  Based  upon  these  two  representa- 
tive problems,  it  would  appear  that  one  oould  expeot  the  application  of 
algorithm  3 to  yield  eolutione  uithin  1 percent  accuracy  in  moot  realistic 
situations.  In  general,  logistics  systems  will  have  more  than  16  SRU's  and 
this  would  further  reduce  the  likelihood  of  a significant  impact  on  the 
solution  arising  from  the  non-convexities  in  the  individual  item  functions 


In  addition  to  algorithms  1,  2,  and  3,  Table  II  has  algorithms  4 and 
5 to  compare  the  E(N0RS)c  values  to  the  approach  when  allocations  are  se- 
lected on  the  basis  of  the  item  that  has  a maximum  for  the  quantity  Z 0 . ./u. . 

i 1 

Algorithm  4 -Starting  at  the  given  input  Lagrange  multiplier  which  is  an 
upper  bound  on  the  constraint,  X^,  and  using  the  item  functions  Uin), 
the  next  allocation  is  selected  by  taking  the  item  that  has  a maxi- 
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bub  for  the  quantity  Y an<i  continuin8  this  Process  until 

there  is  no  item  allocation  reaaining  that  will  not  violate  a constraint. 


Algorithm  5 - The  saae  procedure  as  algoritha  4 except  that  the  starting 
point  is  Ak,  where  k ■ Bax (a- 2,  1). 

The  results  shown  in  Table  II  indicate  that  the  process  of  minimizing  the 
sum  of  expected  backorders  does  not  necessarily  minimize  E(NORS)c.  At 
various  constraint  points  algorithms  4 and/or  5 have  a value  for  E(NORS)c 
significantly  lower  (up  to  31  percent)  than  algorithm  3.  For  the  problem 
that  is  not  presented  in  this  paper  that  has  many  u^  > 1,  the  results  are 
analogous  to  the  exaaple  presented  here.  Also  shown  in  Table  II  are  the 
total  number  of  units  of  stock  allocated  and  the  difference  between  the  solu- 
tion LCC  and  the  problem  constraint  for  LCC. 


Figure  3 shows  a plot  of  the  example  problem  results.  A continuous 

dashed  line  represents  the  expected  backorder  curve  until  it  essentially 

coincides  with  the  E(NORS)  curve.  Thus,  for  LCC  expenditures  beyond 

nc 

$98  M,  the  value  for  E(NORS)nc  could  be  approximated  by  the  expected  back 
order  sum  with  very  little  error,  i.e.,  the  application  of  equation  (14). 
The  bottom  (g)  curve  in  Figure  3 represents  an  approximation  to  the  min- 
imum value  for  E(NORS)c  and  is  selected  as  the  minimum  value  between  the 
solutions  for  algorithms  4 and  5. 


34 


Notice  that  the  g- curve  is  not  convex;  this  is  due  to  either  the  solu- 
tion being  obtained  from  the  non-convex  functions  {5!  } or  the  heuristic 
nature  of  algorithms  4 and  S,  or  both. 

If  one  addresses  the  difficult  problem  of  rigorously  defining  and 
solving  the  problem  of  considering  "some"  cannibalization,  he  soon  finds 
it  mathematically  intractable.  Approaching  the  problem  in  a simplistic  man- 
ner, let  us  define  "the  expected  NORS  where  there  is  x proportion  canniba- 
lization" to  be  a linear  interpolation  between  E(NORS)  and  E(NORS)  ; thus, 

HC  c 

ECNORSV  - xEfNORS)  + (l-x)E(NORS)  . (15) 

x c nc 

Curves  c,  d,  and  e in  Figure  3 show  25,  50  and  75  percent  cannibalization, 
respectively.  The  rationale  supporting  this  approach  to  the  problem  is 
that  the  maintenance  crews  at  the  bases  would  be  instructed  to  consolidate 
the  backorders  from  x proportion  of  the  NORS  vehicles  onto  the  remaining  NORS 
at  the  bases  and  there  is  no  restriction  as  to  which  items  may  be  consolidated. 
Omitted  from  the  problem  is  the  corresponding  increase  in  resupply  time  and 
cost  arising  from  the  consolidation  activity.  This  first  approximation  to  the 
partial  cannibalization  problem  provides  planning  factors  of  sufficient  accuracy 
for  many  applications. 

MODEL  USAGE 

The  principal  logistics  planning  tool  that  evolves  out  of  the  use  of 
the  model  is  the  expected  NORS  versus  LCC  functions  depicted  in  an  illustra- 
tion like  Figure  3.  In  using  the  figure,  the  analyst  wants  to  select  a design 
point  on  the  curve  that  at  least  meets  the  weapon  system's  specification  for 
logistics  availability  and  also  takes  into  consideration  the  marginal  return 
from  additional  logistics  investment,  i.e.,  the  marginal  expenditure  for 
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53  LCC.  to  reduce  the  expected  NORS  by  one  does  not  exceed  the  life  cycle 
i 

cost  to  procure  an  additional  vehicle. 

The  ARIES  solution  at  the  selected  LCC  constraint  defines  the  expected 
backorder  sum  for  each  item.  Target  values  for  the  expected  backorder  sum 
for  each  item  are  a necessary  input  to  the  MS  CAM  model  to  perform  refined 
reliability  and  maintainability  trade-offs  on  each  item.  The  expected  back- 
orders values  for  the  sixteen  items  in  the  example  ranged  from  0.042  to  1.476 
at  a cost  constraint  of  $106M. 

A limitation  with  both  ARIES  and  MSCAM  is  that  items  are  treated  independent- 
ly and,  hence,  the  models  do  not  have  explicit  logic  to  take  into  account  the 
fact  that  several  items  may  use  common  AGE.  Three  approaches  are  identified 
below  to  assist  in  remedying  this  situation. 

Grouping  item  - It  may  be  possible  to  group  several  items  together 
that  use  the  same  AGE  and  have  the  same  maintenance  strategy  and 
other  important  characteristics  and,  thus,  to  treat  these  items  in 
the  model  as  if  they  were  one  item. 

Prorating  AGE  expenses  -It  may  be  possible  to  estimate  the  proportion- 
ate share  of  comson  AGE  usage  by  several  items.  The  input  data  to 
the  model  would  reflect  the  pro  rata  AGE  expenses. 

Attributing  entire  AGE  expense  to  one  item  -It  may  be  appropriate  to 
attribute  the  entire  AGE  expense  to  only  one  item  from  among  several 
that  potentially  could  use  the  same  common  AGE.  This  situation  may 
be  justifiable  if  one  item  has  been  planned  for  AGE-assisted  repair 
in  any  eventuality. 


It  nay  be  necessary  to  iterate  through  the  solution  several  times  to 
arrive  at  the  proper  procedure  for  accomodating  the  comon  AGE  expenses. 

Obviously,  the  types  of  trade-offs  that  can  be  made  with  any  model 
are  variations  of  the  many  combinations  of  input  variables.  For  the  ARIES 
model,  it  is  convenient  to  summarize  the  principal  parameters  for  trade-off 
analyses  according  to  the  following  categorization: 

Reliability 

e Failure  rates 

e Operating  hours  per  base  per  month 

Maintainabi li ty 

e Base  and  depot  AGE 
e Base  and  depot  repair  capability 

1.  Manhours  per  repair 

2.  : Repair  parts  cost  per  repair 

f 

3. /  Repair  turnaround  time 

Availatd  lity 

e System  deployment  concept 
e Expected  item  backorders  across  all  bases 
e Expected  NORS  vehicles 

Item  Charaotevietioe 

e End-item  cost 
e NRTS  and  condeanation  rates 

The  model  has  an  option  to  co^mte  the  unscheduled  maintenance  manpower 
requirements  for  each  item  at  each  base  within  each  skill  category.  Also,  the 


38 


f 


depot  manpower  requirements  and  depot  AGE  utilization  rates  are  calculated 
for  each  item.  Subsequently,  total  base  and  depot  manpower  requirements 
for  all  items  are  summarized  for  convenient  usage  by  the  logistics  planner. 


The  computer  program  for  ARIES  is  written  in  FORTRAN  IV  for  both  the 
CDC  6600  and  the  IBM  370/16S  computers.  The  program  has  three  basic  execu- 
tion sections  that  correspond  to  the  three  previously-discussed  logical 
steps.  The  first  execution  section  processes  the  item  data  to  generate  the 
individual  item  functions  Uin>  • The  common  item  data,  the  item  data, 
and  the  functions  are  stored  on  a permanent  disc  file  for  later  usage. 

The  program  has  a data  management  capability  to  add,  subtract,  or  modify 
items  on  the  permanent  file  to  define  the  system  currently  being  analyzed. 

The  second  execution  section  solves  the  resource  allocation  problem 
subject  to  a single  constraint  by  first  testing  to  see  between  which  two 
input  LagTange  multiplier  constants  the  solution  would  lie.  The  combination 
Lagrange  multiplier  and  final  marginal  allocation  algorithm  reduces  the 
number  of  calculations  over  that  of  simply  using  marginal  allocation  the 
entire  way  for  each  constraint.  The  third  execution  section  solves  for 
the  expected  NORS  and  has  the  option  to  summarize  the  solution  for  each  base 
in  the  system.  A more  detailed  discussion  about  the  ARIES  computer  program 
is  contained  in  reference  11. 

CONCLUSION 


The  ILS  effort  can  be  assisted  in  its  function  to  identify  and  reduce 
I the  logistics  life  cycle  cost  for  a weapon  system  by  employing  an  aggregate 

■ item  model  such  as  described  in  this  paper.  The  model  provides  a general 

1 analytic  treatment  for  optimizing  the  multi-echelon  logistics  system  for 

I 
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unscheduled  maintenance  by  incorporating  the  following  features: 

e Treats  four  alternative  maintenance  strategies  for  an  item 
as  a variable  in  the  problem; 

e Calculates  resource  expenditures  in  terms  of  the  life  cycle  cost; 
e Provides  an  allocation  algorithm  that  produces  solutions  within 
1 percent  accuracy  for  most  realistic  problems; 
e Creates  a permanent  data  base  that  permits  analyses  at  various 

cost  constraints  without  having  to  recalculate  the  entire  data  base; 
e Relates  to  an  operational  measure  of  effectiveness  for  availability- 
expected  NORS; 

e Provides  continuity  with  the  individual  SRU  reliability/maintaina- 
bility trade-off  data  and  model; 
e Computes  base  and  depot  manpower  requirements;  and 
e Provides  useful  resource  allocation  and  performance  data  for  each 
base  and  by  item. 


APPENDIX 


A convenient  recursive  coaputationel  technique  to  reduce  calculations, 
and  to  prevent  underflow  problems  in  the  computer  program  when  either  the 
item  demand  (XT)  or  the  spare  stock  (s)  is  large,  for  the  expected  back- 
order function  is  given  by  the  following  theorem: 


THEOREM.  Let  8 be  the  spare  stock  for  an  item  where  demands  are  described 
by  the  Poisson  distribution  with  mean  arrival  rate  X and  the  resupply  (repair) 
time  is  drawn  from  an  arbitrary  distribution  V(t)  with  mean  T,  then  the 
expression  for  the  recursive  relationship  for  the  expected  number  of  back- 
orders  is 

0(s*l)  - 0(s)-o(s)*y(s), 

where  (16) 

0(0)  - XT,  o(0)  - 1,  y(0)  - e‘XT, 
a(s)  « o(s-l)  -y(s-1)  , and 


Y(s)  ■ y(s-1)XT/s. 

Proof.  Sherbrooke [8]  has  shown  that  0 (s)  can  be  described  by  the  expres- 
sion given  by  equation  (17) 


8(»)  • 


k*» 

£ (k-s)p(k|XT)  , 

k«s*l 


(17) 


where  p(k|XT)  is  the  compound  Poisson  probability  function  for  demands  occur- 
ring during  a mean  resupply  interval,  T.  Using  the  notation  p(k)  for  p(k|XT) 
and  noting  that  for  the  simple  Poisson  p(k)  ■ (XT)ke-XT/k! , then  the  expansion 
of  equation  (17)  yields 


B(s) 


POO  (18) 


k—>  k-s  k=»  k«s 

L kP(k)  - V kp(k)  - s 4 p(k)  * s Y 
k-0  x-0  k=0  k=0 


k-s 

XT  - 4 kpCk) 

k=0 


k-s 

s ♦ s y pck) 
k-0 


k=s-l 

XT  - s ♦ Y (s-k)p(k) 

k-0 

-XT  vk=s_1  k 

4 (s-k)CXT)K/k:  . 

k-0 


- XT  - s ♦ e 


-XT 


Thus,  0(O)-XT  and  6(l)*XT-H-e‘XT-B(0)-a(0)+Y(0)  .where  a(0)  = l and  y(0)-e 
Equation  (19)  results  from  expanding  equation  (18)  with  s+1  being  substituted 

for  s 


6 (s*l)-XT- (s*l)*e'XT  Y kaS (s+l-k) (XT)k/k! 

k-0 


(19) 


-XT-s+e 


-XT  ^*k«s-l 

^k-0 


-XT  _k*s-l, 


(s-k)  (XT)  /k! - (1-e  V*  3 *(XT)VkI] 

i n 


k-0 


♦e~XT(XT)S/s! 


Define  o(s)=o(s-1)-y(s-1)  and  Y(s)-Y(s-1)  XT/s,  then 

B(s+l)=B(s) -a(s)  + Y(s) 

The  terms  a(s)  and  Y(s)  are  readily  updated  as  the  value  of  s is  incre- 
mented by  one.  At  each  value  of  s,  the  quantity  Y(s-l)XT/s  is  tested  against 
the  minimum  permissible  positive  number  for  the  computer, e.  A value  below 
e would  cause  "underflow";  thus,  in  the  computer  program 
y(s)-(y(s-«XT/s  if  Y(s)XT/s>e; 

*>  ’ ^ (0  V v*  otherwise. 


Setting  y(s)  to  zero  when  underflow  would  have  occurred  does  not  affect  the 


t 

value  of  B(s*l)  for  most  realistic  situations  since  c < 10 
scientific  computer  systems 


-75 


in  typical 
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COMPUTER  PROGRAM 
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This  section  presents  a technical  discussion  of  the  computer  program 
for  the  ARIES  model.  The  topics  discussed  include:  program  logic  organization, 
input  data  requirements,  output  displays,  typical  program  usage,  definition 
of  variables,  and  computation  techniques.  The  discussion  centers  around  the 
version  of  the  program  that  inns  on  the  CDC  6600  computer  at  ASD.  There  is 
another,  similar  version  of  the  ARIES  program  that  is  programmed  for  the  IBM 
370/165  computer. 

PROGRAM  LOGIC  ORGANIZATION 

ARIES  is  a large  optimization  system  that  is  composed  of  submodels 
for  performing  specific  computations.  Therefore,  to  describe  the  analysis 
structure  of  the  system,  reference  is  made  to  Figure  4,  which  is  an  overview 
flow  diagram  of  the  system.  The  following  discussion  describes  the  input, 
logic,  and  output  for  each  module. 

Module  1 - MSCAM2 

The  Calspan  MSCAM  model  has  been  modified  (and,  hereafter,  denoted 
MSCAM2)  to  calculate  for  each  item  i,  for  each  of  the  four  alternative  main- 
tenance strategies,  an  appropriate  range  of  expected  backorders  (BO)  and  the 
corresponding  life  cycle  costs,  LCC^  The  input  data  to  Module  1 are  essen- 
tially the  same  as  the  item  input  data  to  the  MSCAM  model.  For  each  of  the 
four  maintenance  strategies,  there  are  decreasing  values  for  the  expected 
number  of  backorders  as  the  LCC^  increases.  The  corresponding  initial  base- 
depot  stockage  levels  are  also  determined  as  a result  of  the  amputation  in 
Module  1. 

Module  2 - Optimum  Item  BO-LCC  Function 

The  inputs  to  Module  2 are  the  data  points  for  the  item  BO-LCC  functions 
for  each  of  the  four  maintenance  strategies.  Module  2 collapses  the  four 
discrete  BO-LCC  functions  to  one  optiaui  function  for  the  item,  foj.  The 

45 


■ 


resultant  function  defines  the  minimum  cost  (also,  the  aalntenance  strategy 
and  the  base-depot  stock  levels)  corresponding  to  a given  value  for  the  item 
expected  backorder  sum. 

After  the  {Cin}  function  is  computed,  the  ratio  6-  from  equation  (6) 
is  calculated  for  each  (BO,  LCC.)  data  point  beyond  the  first.  For  those 
points  that  are  not  on  the  lower  boundary  of  the  convex  hull,  the  value  of  6^ 
is  multiplied  by  -1  to  convey  a non-convexity  signal  and  the  0;  value  within 
one  computer  word.  Using  only  those  points  on  the  lower  convex  hull  boundary, 
the  values  of  are  computed  using  equation  (5).  The  parameters  and  4>^ 
are  used  in  Module  3 with  the  Lagrange  multiplier  and  marginal  allocation 
techniques.  The  output  from  Module  2 and  the  item  input  data  are  stored  on 
random  access  disc  for  later  usage. 

Module  3 - Determination  of  the  Region  Where  a Constraint  ie  Met 


Utilizing  the  principle  of  Lagrange  multipliers,  the  smallest  for 
item  i that  is  equal  to  or  greater  than  specified  Lagrange  multipliers  k^  , 

(k  ■ 1,2,..., IS)  is  noted  and  stored  (with  the  corresponding  S0^  and 
iCci&  f°r  subsequent  use  in  determining  the  total  system  expected  backorders 
and  LCC.  For  each  , the  total  system  expected  backorder  sun  and  LCC  are 
computed  by  equations  (20)  and  (21),  respectively. 


Total  system  — Y SO 

expected  backorders  A*  *A 


Total  system 
LCC 


(20) 

(21) 


Module  3 tests  to  determine  between  which  two  A*  a constraint  is  met 
for  either  total  LCC  or  total  expected  backorders,  or  both;  call  the  two 
bounding  Lagrange  multipliers  kmf  and  km . Module  3 prints  out  the  total 
LCC  and  the  total  system  expected  backorders  for  each  Aj  . 
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Module  4 - Compute  Solution  at  Constraint 

Once  the  region  has  been  defined  where  either  a cost  or  backorder 
constraint  is  met,  the  technique  of  marginal  allocation  is  employed  to  reach 
the  solution  at  the  constraint.  The  following  procedure  is  performed: 

(1)  Initialize  the  sum  of  system  backorders  and  total  cost 

by  summing  the  individual  item  backorders  and  cost  at 
where  p • , and  \m  is  where  a constraint  is 

met  in  Nodule  3.* 

(2)  Using  marginal  allocation,  select  the  item  which  has  the 
largest  ratio  of  marginal  decrease  in  expected  backorders 
divided  by  the  marginal  increase  in  LCC  to  obtain  an  expected 
backorder  reduction  (i.e.,0/  from  equation  (5)). 

(3)  Test  to  see  if  a constraint  would  be  violated:  if  yes, 

go  to  step  (5);  otherwise,  go  to  step  (4). 

(4)  Update  the  total  system  backorders  and  cost  by  adding 
the  additional  unit  of  the  item  selected  in  step  (2) 

to  the  system.  If  using  , go  to  step  (2);  otherwise, 
go  to  step  (5). 

(5)  Using  marginal  allocation,  select  the  item  which  has 
the  largest  value  (i.e.,  from  equation  (6)). 

(6)  Repeat  steps  (3)  through  (S)  until  there  are  no  more 
feasible  item  allocations  that  can  be  added  to  the 
system  with  a constraint  still  not  violated. 


The  program  actually  has  five  different  variations  of  the  marginal  allocation 
technique.  (Refer  to  Section  IV  for  the  definitions  and  to  Card  Group  IX 
(Table  IV)  for  the  input  data  format.) 
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The  benefit  of  using  the  Lagrange  Multiplier  technique  in  conjunction 
with  Marginal  allocation  is  that  a large  reduction  can  be  made  in  the  nusber 
of  calculations  while  a solution  can  be  approached  that  is  very  close  to 
optimal  for  the  constraint. 

Module  5 - Calculate  Base  and  Item-Depot  Swmariee 

The  solution  computed  in  Nodule  4 defines  the  item  maintenance 
strategies  and  the  base  and  depot  stock  levels  [s-  ; for  all  < and 

j J . Using  the  item  input  data  and  the  item  problem  solution  data,  modified 
versions  of  MSCAM-type  submodels  are  employed  to  compute  useful  summary  data 
for  each  base  and  also  at  the  depot.  Tables  VI  and  VII  (which  appear  later 
in  this  section)  illustrate  how  the  output  data  would  be  arranged  for  the  base 
and  item-depot  susnaries,  respectively.  The  expected  number  of  NORS  missiles 
is  computed  for  the  case  when  complete  cannibalization  is  assumed  and  for  the 
case  when  there  is  no  cannibalization. 

Subroutine  Structure 

ARIES  is  composed  of  17  subroutines.  Figure  5 shows  how  the  sub- 
routines interface.  Table  III  sumarizes  the  function  of  each  subroutine  and 
gives  its  location  in  the  program— i.e.,  overlay,  overview  logic  module,  and 
program  execution  section  (discussed  in  the  next  subsection  below). 

Program  Execution  Sections  and  Disc  Data  Piles 

A convenient  way  to  visualize  the  ARIES  computer  program  is  to 
consider  the  program  as  consisting  of  three  execution  sections.  The  first 
section  processes  the  item  data  to  generate  the  individual  expected  backorder- 
LCC  functions,  . Section  I corresponds  to  principal  subroutines  STOCK 

and  CAL.  The  item  data  are  stored  on  disc  file  8 and  can  be  treated  as  a 

permanent  data  file. 
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I 


Computer  Program  Subroutines 


A 


Execution  Section  II  corresponds  to  subroutine  SETA  in  the  program. 
SETA  reads  file  8 twice.  SETA  generates  two  temporary  disc  files  for  inter- 
ned! ate  calculations.  File  9 is  used  to  store  the  data  from  file  8 for  the 
region  where  a constraint  is  met.  File  10  is  used  to  store  the  table  of 
index  values  that  define  the  individual  item  solutions  for  the  fifteen 
Lagrange  multiplier  values. 


Execution  Section  III  corresponds  to  subroutines  NORS,  BOB,  and  SPOT. 
NORS  reads  file  9 once.  NORS  generates  temporary  file  11.  File  11  is  used 
to  construct  the  output  tables. 

INPUT  DATA 

Table  IV  defines  the  ARIES  input  data  requirements.  The  first  data 
card  defines  the  value  of  EXMODE  and,  hence,  which  sections  of  the  program  are 
to  be  executed. 


EXMODE 


SECTIONS  EXECUTED 


I only 

I and  II 

II  and  III 
II  only 

II  and  III 


Whenever  Section  I is  executed,  a file  8 is  either  created  (if  there 

is  no  previous  file  8)  or  extended  with  the  data  for  new  items  that  are  under 

calculation  with  the  current  computer  run.  When  there  is  no  pre-existing  file 
8,  set  I0M0DE  to  0;  when  there  is  an  existing  file  8,  set  IOMODE  to  1.  If 

EXMODE  equals  4,  5,  or  6,  a pre-existing  file  8 is  assumed. 

EXMODE >6  differs  from  4 by  the  fact  that,  with  6,  the  program  assumes 
that  a constraint  problem  has  been  solved  prior  to  the  case  with  EXMODE-6  (i.e., 
during  the  same  computer  run),  with  the  same  Lagrange  multipliers,  the  same 


zx.ii  mmmii 


TABLE  OF  (Cont.) 

INPUT  DATA  FORMAT  FOR  THE  ARIES  MODEL 


J 

1 

1 


Ctrd  Group  III  ICont.l 

T n - Card  W malted) 


OataPM 

far  at 

Variable 

Brief  Explanation 

1-10 

F10.J 

OOTO.  J) 

lama  a*  above 

11-JO 

F10.J 

OI0T(l.  J> 

Jama  or  above 

21-30 

F10.3 

OSTO,  J) 

Jama  at  above 

31-40 

F10.3 

DilTtl.  J) 

lame  at  above 

CM  Group  IV  (2  cMs): 

Fermat  FMT3 

For  at 

Variable 

Brief  Explanation 

FMT3 

FOAC(I) 

For  dapot  only  and  no  repair  ttratapim.  the  fraction 
el  bate  AGE  aoat  incurred  for  equipment  type  1; 
107 

FMT3 

FRF 

Failure  rate  factor  to  multiple  the  input  failure 
rata  XX(21> 

FMT3 

FWF 

10- Year  uniform  tariot  prevent  worth  factor 

CM  Group  V (l  card):  FOTO 


ted  Oreupt  VI  dw^i  VIII  me  rap mm*  ft 

CM  Group  VI  (UmmMy  3 cm*):  Fermat 


tea  repair  aatt.  0 par  man. 
Dapot  repair  eoat.  $ par  mai 
Mpptap  aoat,  l/*on /mile 


Variabla 

XX(B)  to  XXI7) 
XX  (01 

Ham  10 
Factor  denol 

Brief  Explanation 
inp  the  number  of  timet  the  bam 

xxm 

repair  mat 
to  dataat/ 

t heart  axoaad  the  bate  ^man-bour^  ^ 
r 1 to  J in  ted  Group  III 

XXOO) 

Number  af  a 

nd  team  uniti  par  vehicle 

xxnii 

AGE  inaaaOa 

don  and  tet-up  coat,  0 

XXI12) 

Number  el  | 

npm  ad  taahniaal  data 

XXI13I 

Annual  open 

•dap  and  malxtananae  aoat  par  AGE. 

XX(14I 

Coat  af  an  a 

nd-HMt,  0 

XXflEI 

End  Ham  tb 

Ippinp  aaaipbt  ft 

Bar 


TABLE  DT  (Corn.) 

INPUT  OATA  FORMAT  FOR  THE  ARIES  MODEL 


Coni  Group  VI I Com  / 


XXIIt) 

XXI17) 

xxnai 

XXI19) 

XXI30I 

XXI21) 

XX(22I 

XXI23) 

XX(M) 

XX(tt) 

XX(2ti 

XXI27) 

XX(») 

xxiai 

XXI30) 

XXI31) 


TABLE  IS  (Cont.) 

INPUT  DATA  FORMAT  FOR  THE  ARIES  MODEL 


Cmd  Group  IX  (2  curds/  For  nut  (2E16.7.  416/2014 


Og|yHD 

fount 

Variable 

ff.l  .1  m * -*  — 

cntT  c xpnniiion 

1-16 

016.7 

OMAX 

C«t  constraint 

16-30 

016.7 

on 

Bachonhr  conatraint 

31-36 

16 

MO 

Maximum  numbor  of  data  point*  par  itam 
on  fit*  8;  i 300 

36-40 

16 

OLCTNO 

Number  of  items  ratMinp  on  «•*  6 that  art 
to  ba  ddated  from  tha  currant  execution 
of  SETA;  £.10 

41-46 

16 

iorr 

Allocation  algorithm;  IOFT  - 1,  2.  3.4.  or  6 Saa 
Section  IV  for  definition*.  Generally.  IOPT  - 3 
i*  u*ad. 

46-60 

16 

IZ 

For  IZ-IOPT,  fame  (ummarie*  are  printad  out 
For  IZdlOPT.  bate  Mimmary  printout  to  *uppra*aed 

6mm0«** 

OdImimi 

Format 

Variable 

Brief  Explanation 

1-4 

M 

NV(I) 

Number  of  eohido*  at  ba*a  1 

66 

14 

NV(II 

Number  of  Vahid  a*  at  taaw  2 

0-13 

•4 

NV<1) 

Number  of  vefude*  at  ban  3 

00-66 

14 

Nvin 

Number  of  vehicle*  at  ban  17 

Curd  Group  X (3  curds): 

Format  (60107) 

* 

CM 

fwnt 

VbHbMb 

Brief  Explanation 

1-16 

016.7 

LMO) 

Layanpa  multiplier  value*  in  strictly  monotontcaUy 
daarntinp  order 

10-30 

010.7 

LNM2) 

31-40 

0107 

« 

40-00 

0107 

. 

01-10 

0107 

• 

N OLCTNO  to  Mrs.  4m k 

Card  Orcup  XI  to  ildppad 

Oust  Group  XI  (1  curd  for  tuck  Item  deleted/ 

Format  (3A4) 

ry 

NrfM 

wniPi 

1 

Mi 

1 

1-12 

S M 

(DO  LIT  (1.  J). 
1*1,31 

Ham -10  of  Ham  rnidht|  on  f Ha  0 to  be  dal  and  Horn 
aurrent  execution  of  SETA  and  NOM 

End  of  DATA  DICK 
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number  of  vehicles  per  base,  and  the  same  items  that  are  going  to  be  used  for 
the  current  constraint  problem. 

For  the  CDC  6600  computer,  the  control  card  "CATALOG,  TAPE8,  ARIES, 

CY  ■ RP  ■ 999."  is  used  to  retain  the  original  generation  of  a file  8.  The 
control  card  "EXTEND,  TAPE8."  is  used  to  extend  a pre-existing  file  8 with 
modified  item  data  or  additional  item  data  and  to  retain  it  permanently. 

The  capability  exists  to  suppress  items  on  file  8 from  the  optimiza- 
tion calculations  in  Section  II  and  III  by  use  of  the  input  variables  DLETNO 
and  DELET.  The  parameter  DLETNO  (Card  Group  IX)  indicates  the  number  of 
records  of  file  8 that  are  to  be  deleted;  the  item-IDs  of  the  items  to  be 
deleted  are  placed  in  array  DELET  (Card  Group  XI).  A maximum  of  ten  items 
from  file  8 may  be  deleted  on  a given  computer  run. 

Unlike  the  modifications  which  result  from  extending  file  8,  use  of 
DLETNO  and  DELET  does  not  cause  a "physical"  change  in  file  8 (i.e.,  there 
are  no  records  actually  removed  from  the  file).  The  actual  effect  is  that 
records  with  item-IDs  appearing  in  DELET  are  ignored  during  the  current  calcu- 
lations in  SETA;  they  remain  in  file  8 and  are  accessible  for  subsequent 
executions  of  SETA. 

TYPICAL  PROGRAM  USAGE 

ARIES  is  written  in  extended  FORTRAN  IV  language.  The  program  uses 
106.K  central  memory  words  and  takes  approximately  18  seconds  to  compile. 

o 

The  central  processor  time  to  generate  the  item  expected  backorder- LCC  functions 
for  the  example  problem  discussed  in  Section  IV  is  approximately  90  seconds. 

It  takes  0.062  seconds  of  c.p.  time  to  execute  Sections  II  and  III  of  the 
program  for  a constraint  in  the  example  problem;  the  1/0  time  is  approximately 
31  seconds.  The  1/0  time  is  relatively  high,  since  the  CDC  algorithm  for  disc 
1/0  time  incorporates  a large  fixed  time  for  each  instance  in  which  the  disc 
is  called  in  addition  to  the  variable  time  relating  to  the  number  of  words 
transferred. 
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A typical  way  to  use  ARIES  would  be  to  generate  file  8 with  EXMODE-1, 
a low  backorder  constraint,  and  a cost  constraint  that  would  likely  fall  within 
the  range  of  the  data  for  the  15  Lagrange  multiplier  points.  Upon  examining 
the  output,  note  the  upper  and  lower  limits  of  the  total  system  LCC  for  the 
range  of  15  Lagrange  multipliers.  Using  EXM0DE=4  for  the  first  constraint 
and  EXMODE-6  thereafter,  submit  another  run  with  sets  of  data  for  total  cost 
constraints  spaced  over  the  noted  cost  interval  to  derive  the  curve  for 
expected  NORS  (availability  measurement)  versus  LCC.  The  program  is  written 
to  execute  new  data  sets  until  an  end-of-file  card  is  encountered. 

The  item  input  data  cards  (Card  Groups  V through  VIII)  are  same  as 
those  for  the  Calspan  version  of  the  MSCAM  model  with  the  exception  that  XX (9) 
for  MSCAM  is  the  "expected  sum  of  item  backorders  across  all  bases"  instead 
of  the  "depot  index  corresponding  to  J in  Card  Group  III."  This  high  degree 
of  commonality  facilitates  continuity  of  operation  between  the  two  models. 

DEFINITION  OF  VARIABLES 

Since  ARIES  is  a large  program,  it  would  be  too  voluminous  to  define 
here  each  FORTRAN  variable.  Table  V presents  some  of  the  principal  variable 
definitions,  along  with  their  locations  in  the  program.  The  input  data  require- 
ments (Table  IV)  also  define  program  variables. 

OUTPUT  DISPLAYS 

ARIES  is  written  to  print  out  several  formal  tables  and  several  inter- 
mediate calculations  which  are  helpful  for  checking  that  the  results  are 
logical  and  consistent.  The  content  of  two  of  the  formal  tables  is  illustrated 
by  Tables  VI  and  VII,  which  show  the  base  and  item-depot  summaries,  respec- 
tively. 

In  each  table,  the  term  "Not  Operationally  Ready  - Supply"  (NORS)  is 
listed.  For  the  base  summary  output  table,  the  expected  number  of  NORS 
missiles  attributable  to  a particular  item  is  defined  to  be  the  expected  back- 
order value,  divided  by  the  number  of  units  of  that  item  which  are  used  on  a 
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I 
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1 

I 

I 

1 

1 
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TABLE  T 

DEFINITIONS  OF  SOME  OF  THE  VARIABLES  APPEARING  IN  ARIES 


FORTRAN  NAME 


DEFINITION 


EXECUTION  SECTION 


NITEMS 


NDPTS 


PROG(I) 

OST(l,  J) 

DIST(I,  J) 
BRCST.XX(I) 
DRCST,  XX(2) 
SHPCST,  XX 13) 
PAKCST,  XX(4) 
SETCST,  XX(11) 
DATCST,  XXII 2) 
OPCST,  XXI13) 

ENDCST.  XXII 4) 
ENDWT,  XXI15) 
EXNRTS,  XXII 6) 
CNDMX,  XXII 7) 
BRUME,  XXII 8) 
BMANHR,  XX|1 9) 
DMANHR,  XXI20) 
FAILRT,  XXI21I 
BAGEWT,  XXI 22) 
DAQCST.  XXI23) 
DLOTME,  XXI24) 
BPTCST,  XXI26) 
DPTCST,  XXI26) 
NSKIL,  XXI27) 

BAQCST,  XXI28) 
DRTIME,  XXI 29) 
DAQEHR,  XXI30) 
NITS 


Execution  option  parameter 

Numbar  of  Rants 

Numbar  of  basat 

Numbar  of  dapots 

Input  device  numbar 

Control  to  stop  calculations  in  STOCK 

Control  to  stop  calculations  in  STOCK 

Disk  input  option  paramatar 

Total  basa  oparating  hours  par  month 

Ordar  and  shipping  time 

Distanca  from  basa  I to  dspot  J 

Basa  rapair  cost,  S par  man-hour 

Dspot  rapair  cost,  S par  man-hour 

Shipping  eost/ton/mile,  t 

Packaging  oost/shipmant.  $ 

Installation  and  AGE  aatup  cost 
Numbar  of  pagas  of  technical  data 

Annual  operating  and  maintenance  cost 
par  unit  AGE 

Cost  of  end-Ram 
Shipping  weight  of  end-Ram 
Expected  NRTS  rata 


Average  basa  rapair  time,  days 
Man-hours  par  )ob  to  fix  at  bate 
Man-hours  par  Job  to  fix  at  depot 
Failures  par  100  effective  operating  hours 
Weight  of  bate- AGE 
Cost  of  depot- AGE 
Procurement  lead  time 
Rapair  parts  ooet  at  bate 
Repair  parts  oott  at  depot 

o.a_  . a_  _ M.  -a-sat-  a— s_  « 

FvunHJvr  pt  dfnirvfn  wiiii  ip  rwpmr  mil 

Ram  at  basa 

Cost  of  bate- AGE 

Average  repair  time  at  depot 

Average  hr /repair  Depot  AGE  la  required 


Index  no.  for  ARCAL  do  loop 


EXECUTION  SECTION 


TABLE  V (coot.) 

DEFINITIONS  OF  SOME  OF  THE  VARIABLES  APPEARING  IN  ARIES 


CODE  A(1 ),  CODE  BO) 

PNTSKL 

DEMLT 

BOFM).  J - 1.  2 

BOOM) 

MP 

M 

DISTO,  J) 

GSM.  1) 

GS(J,  2) 


NBO,  U 

TABLE 

BT(K) 

LCC(K) 


Percentage  repair  par  labor  skid 
Avaraga  demand  during  lead  time 

Backorders  at  random  point  in  time;  BOF(I)  it  back- 
orders with  one  lam  unit  of  stock  than  BOF(2) 

Average  delay  at  depot  given  depot  stock  U-1) 

Upper  limit  for  item  total  stock 

MP-1 

Distance  (miles)  between  depot  J and  base  I 
Sum  of  base  level  stock  for  cate  J ■ 1, ....  4 
Depot  stock  level  for  cere  J*  I, ....  4 
Backorder  array  as  returned  by  STOCK 
Cost  array  as  returned  by  STOCK 
Backorder  array  as  returned  by  CAL 
Coat  array  as  returned  by  CAL 
Backorder  array  passed  from  SETA  to  NORS 
Cost  array  pasted  from  SETA  to  NORS 
fUrkorrir  imv  at  rtturnBd  from  BOB 

e^^^^werv  s^ev  w e o y e ewm  v^^as  a a tress  ^^tesr 

Backorder  array  at  base  I at  returned  from  SPOT 


EN 

XMNYR(K) 


<t>  array 

Index  counter  for  file  8 

Index  counter  for  file  8 record  number 

Current  value  of  minimum  coat  that  satisfies 
and  *2 

Diagonal  array  of  backorders  from  the  base  da 
allocation  matrix 

Number  of  units  of  stock  at  bate  I for  L-1 
total  units  of  stock 

Array  with  ID-namet  of  the  Hems  on  file  8 

Total  expected  backorders  at  Lagrange 
multiplier  K 

Total  LCC  at  Lagrange  multiplier  K 

Total  expected  backorders  at  constraint  point 

Total  LCC  at  constraint  point 

Index  counter  denoting  the  position  on  file  8 
whore  a Lagrange  multiplier  point  is  mat 
for  the  item  under  current  calculation 

Expected  NORS  assuming  oannibaliaation 

Expected  NORS  assuming  no  cennibeiisation 

Number  of  man-years  for  labor  skill  K 


missile.  The  number  of  NORS  missiles  for  each  base  is  equal  to  the  maximum 
value  for  an  "item  NORS"  when  parts  are  consolidated  through  the  practice  of 
cannibalization.  NORS  missiles  will  be  constrained  by  the  item  in  shortest 
supply.  The  expected  NORS  index  for  each  item  listed  in  the  item-depot 
summary  is  the  item  expected  backorder  value,  divided  by  the  number  of  units 
of  the  item  that  are  used  on  the  missile.  The  expected  NORS  value  listed  at 
the  bottom  of  the  table  is  the  simi  of  the  maximum  expected  NORS  values  from 
all  bases  (the  sum  of  the  "total"  NORS  from  the  base  summaries). 


The  value  for  expected  NORS  missiles  assuming  no  cannibalization  is 
also  computed  and  printed  out  separately. 


A table  is  printed  out  to  define  the  index  counter  value  for  the 
BO-LCC  function  for  the  solution  for  each  item  at  each  Lagrange  multiplier 
point.  This  table  provides  a useful  indication  of  how  the  item  stockage  levels 
are  distributed  over  the  range  of  system  expenditure  levels. 


A table  is  printed  out  to  show  the  total  system  expected  backorders 
and  the  LCC  for  each  of  the  IS  Lagrange  multiplier  points.  This  table  is 
useful  for  determining  how  total  cost  constraints  should  be  set  to  derive  the 
total  expected  NORS  versus  total  LCC  function. 


COMPUTATION  TECHNIQUES 


This  subsection  discusses  some  of  the  special  computation  techniques 
that  are  employed  in  the  ARIES  program. 


Subroutines  STOCK  and  SPO 


For  each  item,  the  BO-LCC  function  must  be  computed  for  each  of  the 
four  alternative  maintenance  strategies.  SPO  computes  the  diagonal  of  a 
matrix, (akf),  where  each  element  is  an  optimal  allocation  of  k units  of 
stock  among  all  the  bases  and  p units  of  stock  at  the  depot.  The  allocation 
is  performed  by  iteratively  selecting  the  next  unit  of  stock  for  a base  that 


will  maximize  the  reduction  in  the  sum  of  expected  backorders  at  all  the  bases. 
STOCK  searches  the  diagonal  of  constant  total  stock  (i.e.,  k+p  - constant)  for 
the  base-depot  stockage  level  combination  that  will  provide  the  minimum  ex- 
pected backorder  sum.  The  program  terminates  the  calculations  for  an  item 
maintenance  strategy  when  any  of  the  following  conditions  are  satisfied: 

(1)  The  strategy  under  current  calculation  is  the  first 
one  for  the  item  (i.e.,  base-depot),  and  the  expected 
backorder  value  is  a €,  and  s 6Z  , where  €1  and  are 
input  parameters;  or 

(2)  The  total  number  of  units  of  stock  > 300;  or 

(3)  The  strategy  under  calculation  is  not  the  first  one 
(i.e.,  it  is  D,  B,  or  NR),  and  there  is  a previously 
calculated  strategy  that  has  a maximum  LCC  (at  the 
point  of  satisfying  «,  and  ) which  is  lower  than 
the  current  value  for  LCC  for  the  strategy  under 
computation.  Stated  analytically,  there  is  a point 
beyond  which  the  strategy  is  dominated  by  at  least 
one  other  strategy. 

Subroutines  PDIST  and  PD2 

An  efficient  recursive  computational  technique  is  used  in  PDIST  and 
PD2  to  compute  the  expected  backorder  values.  This  procedure  is  discussed 
in  the  appendix  of  Section  IV. 

Subroutine  HULL 

The  basic  procedure  employed  to  derive  the  convexified  item  expected 
backorder-LCC  function,  { , involves  the  iterative  testing  (starting  from  the 
leftmost  point)  to  see  where  the  next  point  to  the  right  lies  that  has  the 
minimus  slope  between  it  and  the  last  point  found  on  the  convex  hull.  If  the 
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point  found  is  not  adjacent,  then  a concave  region  exists;  this  region  is 
identified  in  the  prograa  by  Multiplying  the  in  this  region  by  -1. 

PROGRAM  CHECKOUT 

The  ARIES  computer  program  was  checked  out  with  two  extensive  test 
cases.  The  first  test  case  is  the  example  problem  discussed  in  Section  IV. 
The  input  data  and  sample  output  displays  for  this  test  case  are  given  in 
Supplement  A. 

The  second  test  case  involved  SRUs  that  have  more  than  one  unit 
used  per  missile  (i.e.,  u->f  for  some* ).  The  input  data  and  sample  output 
displays  for  this  test  case  are  presented  in  Supplement  B. 


SuDotamant  A 

EXAMPLE  PROBLEM  INPUT  AND  OUTPUT  DATA 

This  supplement  presents  the  input  data  and  sample  output  displays 
for  the  example  problem  discussed  in  Section  IV. 


INPUT  DATA 


The  input  data  for  the  example  problem  are  summarized  below. 


Card  Group  I 

EXMODE- 1 NDPTS-1  EPS2*2.5xlO'2 

NITEMS-16  NINP-S  IOMODE-O 

NBASE-17  EPS1=5x10-8 


Card  Group  II 


(6F10.2) 

(3A4/(6F10.0) 

(7F10.2) 


Card  Group  III 


I_ 

PROGCI) 

OST(I.l) 

DIST(I.l) 

1 

490. 

25. 

850. 

2 

490. 

25. 

900. 

3 

490. 

25. 

850. 

4 

490. 

25. 

700. 

S 

252. 

25. 

1500. 

6 

2S2. 

25. 

1100. 

7 

252. 

25. 

250. 

8 

252. 

25. 

250. 

9 

490. 

25. 

800. 

10 

490. 

25. 

700. 

11 

252. 

25. 

1200. 

12 

252. 

25. 

750. 

13 

490. 

25. 

1200. 

14 

490. 

25. 

600. 

IS 

252. 

25. 

400. 

16 

252. 

2S. 

900. 

17 

252. 

25. 

1200. 

66 


Card  Group  IV 


I 

I 

I 


PBAC(I) 

FRF»3. 


I 

0.90 

PWF  - 10. 


2 

0.95 


3 

0.10 


4 

0.20 


5 

0.05 


Card  Group  V 


l 

XX(D 


1_ 

11. 


2 

22. 


3 

0.5 


4 

200. 


Card  Group  VI 


See  Table  VIII 


Card  Groups  VII  and  VIII 

100AA  ) £or  each  itejn 

1.0  J 

Card  Group  IX 


DMAX»93.5xl06,  94xl06,...,  llOxlO6;  DMAX  is  varied  for  each  new 

constraint  problea. 


EPS 

■ 1.0,  MP  - 

300; 

DLETNO  » 

0;  IOPT 

- 3;  IZ 

» 0 

l 

Nmi 

1 

NVHI 

_I 

NV(I) 

_I 

Mil 

1 

70 

6 

36 

10 

70 

14 

70 

2 

70 

7 

36 

11 

36 

IS 

36 

3 

70 

8 

36 

12 

36 

16 

36 

4 

70 

9 

70 

13 

70 

17 

36 

S 

36 

67 


: 


! 

5 


I 

I 

I 

} 


5 


m 

S | 

SIB 

DBI 

a 

a 

B 

BI 

pi 

S * 1 
a.,,  a 

■B 

■ 

01 

Q 

HI 

Si 

i 

a 

B 

ill 

S 

K>  ] 

■fl 

SB 

■ 

* i 

Si 

Si 

i 

a 

B 

!S 

11 

mm 

in 

Bn 

H 

* i 

Si 

III 

a 

B 

m 

11 

M 

fuA 

K>- 

H 

1 * 
s 

g a 

e 

1471 

67 

11333 

HI 

i 

a 

*1 

s 

S 

mm 

£ - a 

■ 
v&C*  ' 

« 

i 

8 ! 

5 6 S 

8 R 8 

i 

a 

*1 

ib 

S 

Is" 

C 

8 \ 

m — 

| 

1*8 

I b n 

S 

5 

i 

a 

D 

BB 

B 

1 

I • 
* 1 

id 

i 

8*8 

|M 

5 

5 

i 

s 

fl[ 

s 

B 

i 

g 

i 

BiD 

IBB 

« ! 

a 

?! 

« 

s 

B 

1 

HI 

0 

SD 

SB 

i 

a 

Bi 

II 

B 

1 

i 

SD 

SB 

5 

i 

a 

Ml 

s 

B 

IB 

HI 

B 

SD 

SB 

i 

a 

wm 

■ 

s 

bi: 

IB 

§ 

SD 

SI 

■a 

i 

a 

IS 

■ 

BI 

bi: 

PI 

■■ 

BID 

IS 

5 3 

i 

a 

M 

IB 

IK 

BID 

DBB 

i 

a 

B 

BB 

BE 

i 

BSD 

IBB 

i 

a 

B 

S 

BE 

i 

iff 

5 

8 

5 5 3 

a aS 

1 

t i 
81 

i i ] 

is! 

II 

II 

II! 

11! 

II 

Card  Group  X 


mu 

MU 

_I 

LM(I) 

1 

9.16227766x10"! 

6 

6 . 309573445x10"^ 

11 

1.995262315x10"^ 

2 

1.99526231x10"® 

7 

5. 011872336x10 

12 

1.584893192x10-6 

3 

1. 258925412x10"® 

8 

3.981071706x10-6 

13 

l.OxlO"6 

4 

l.xlO-5 

9 

3. 16227x10"® 

14 

6.630957x10-7 

5 

7. 943282347x10”^ 

10 

2.51186xl0"6 

15 

3.162278xl0-7 

Card  Group  XI 

Not  needed  for  this  problem. 

SAMPLE  OUTPUT  DISPLAYS 

Table  IX,  X,  and  XI  show  sample  printout  displays  for  the  example 
problem.  Table  IX  shows  the  total  system  sum  of  expected  backorders  and  LCC 
corresponding  to  the  15  Lagrange  multipliers. 

TABLE  XC 

TOTAL  SYSTEM  SUM  OF  EXPECTED  BACKORDERS  AND  LCC 
FOR  THE  EXAMPLE  PROBLEM 


K 

1 

2 


LAM90MK) 

• 8l62278E*M 
•19952S2F-G4 
• 1258925E-34 
•l000003E-:4 
•7943282E-G5 
•6309573E-3S 
•5011872E-05 
•3981072E-05 
•31S2278E-35 
•2511886F-05 
•1995262F-05 
•1584893E-0* 

• 1JG0  3Q0  F-35 
•6309573E *36 
•3162278E-06 


3TKI 

• 1S26232EO? 
.A171959E02 
,«5741^4E»32 
•4319U*  8E*32 

• 4334299E02 

• 251884CE02 
. 2 J45734E0  2 
.153236 9E*J2 

• 1338587Ef’j2 

• 107423  7£f 92 

.ajraiaaeaJt 

• 513432  6C+  3 1 
•3499816E*0l 

• 276905 3 E+01 
.8754266EO0 


L22K» 

• 9349321E4-9A 
•9517S73E*08 
•9665l53Et08 

• 973355SE*9  8 

• 982C220EO8 

• 103  35  B1E>0  9 
•1106497EfL9 
•1025263E»39 
.1030477E>09 
•104C251E+Q9 

• 1052563EO9 
•106B334E*09 
•1079692F*C9 
•1088624E*99 
.11297  35E+09 


Table  X is  a tabulation  of  the  KJ  variable  for  each  item  for  each 
of  the  15  Lagrange  multipliers.  KJ  is  an  index  counter  corresponding  to  tti 
for  item  i in  Equation  (6).  Thus,  KJ  reflects  the  sequence  number  for  the 
original  item  functions  £ that  represent  the  item  solution  at  each 
Lagrange  multiplier  value.  The  value  of  KJ  is  often  close  to  representing 
to  total  number  of  units  of  stock  for  the  item  that  is  allocated  between  the 
bases  and  the  depot;  however,  it  may  be  different  when  the  item  maintenance 
strategy  changes  as  the  item  resource  expenditure  level  is  increased. 

Table  XI  shows  the  item  and  depot (s)  summary  printout  for  a LCC 
constraint  of  $98  M. 
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Supplement  B 


TEST  CASE  WITH  SOME  SRUt  THAT  REQUIRE  MORE  THAN 
ONE  UNIT  PER  VEHICLE 

This  supplement  presents  the  input  data,  solution  results  summary, 
and  sample  output  displays  for  a model  test  case  that  has  some  SRUs  that 
require  more  than  one  unit  per  vehicle. 

INPUT  DATA 


The  input  data  for  this  test  case  are  summarized  below. 
Card  Group  I 


EXMODE=l 

NITEMS=15 

NBASE=17 


NDPTS=3  ESP2=0.025 

NINP=1  IMODE=0 

ESP1=5x10-6 


Card  Group  II 


(5F10.2) 

(3A4/F12.0,  2F4.0,  6F10.2/7F10.2/8F10.2) 
(7F10.2) 


Card  Group  III 
I PROG(I) 

OST(I.l) 

DISTCI.l) 

0ST(I,2] 

DISTCI.2) 

0ST(I,3) 

DISTfl 

1 

160 

25 

850 

30 

900 

20 

800 

2 

160 

25 

900 

30 

1000 

20 

3000 

3 

160 

25 

850 

30 

1000 

20 

450 

4 

160 

25 

700 

30 

800 

20 

700 

S 

160 

25 

1500 

30 

700 

20 

800 

• 

160 

25 

1100 

30 

600 

20 

900 

7 

160 

2S 

250 

30 

2000 

20 

700 

• 

160 

25 

2S0 

30 

1000 

20 

800 

• 

160 

2S 

•00 

30 

1500 

20 

900 

M 

160 

2$ 

750 

30 

600 

20 

700 

M 

MO 

n 

1200 

30 

800 

20 

1200 

If 

MO 

n 

710 

30 

700 

20 

1000 

n 

IJM 

30 

1000 

20 

800 

• 

m 

600 

» 

soo 

20 

700 

4 

m 

400 

m 

1000 

20 

900 

m 

• 

tm 

m 

1200 

20 

500 

^0 

• 

iMO 

» 

•00 

20 

1100 

*» 


! 

1 

1 

P 

Card  Group  IV 

1 

r 

III 

3 4 

5 

1 

PBAC(I)  0.95  0.5 

0.2  0.8 

0.1 

1 

FRF=2.  ; PWF=10. 

J 

Ti 

Card  Group  V 

] 

III 

3 4 

1 

XX(I)  9.  22. 

.5  100. 

■ 

Card  Group  VI 

See  Table  XII 

1 

i 

•• 

Card  Groups  VII  and  VIII 

] 

100AA  \ £or  item 

1.0  J 

1 

Card  Group  IX 

1 

DMAX-71.X106,...,  $9xl06 

; DMAX  is  varied  for  each  new 
constraint  problem. 

EPS  = 1.0;  MP  = 300;  DLETNO  = 0;  IOPT  = 3; 

IZ  = 

o; 

NV(I)  = 60  for  I ■ 1,  2, 

....  17 

] 

Card  Group  X 

] 

I LM(I)  I 

MI) 

_I 

LMCII 

Tl 

1 9.16227766x10“!  6 

6.309573445x10-6 

11 

1 .995262315x10“^ 

J 

2 1.99526231x10“^  7 

5.011872336x10-6 

12 

1.584893192x10-6 

Tl 

3 1 . 258925412x10"^  8 

3.981071706xl0“6 

13 

l.OxlO-6 

1 

4 l.xlO"5  9 

3.16227x10"^ 

14 

6.630957X10"7 

5 7.943282347x10-6  10 

2.51186x10"^ 

15 

3.162278xl0“7 

] 

74 

] 

# 

] 

J 

rJI 

* m 

- - -1 

jfc. - . 

TABLE  ZD 

TEST  CASE  (WITH  SOME  Uj  > 1)  INPUT  DATA  FOR  THE  XX-ARRAY  (CARD  GROUP  TD) 


Card  Group  XI 


Not  needed  for  this  problem. 

ANALYSIS  RESULTS  SUMMARY 

The  analysis  results  from  solving  a range  of  LCC  constraints  problems 
for  the  test  case  are  summarized  in  Tables  XIII  and  XIV.  The  five  allocation 
algorithms  and  column  headings  in  Tables  XIII  and  XIV  are  the  same  as  those 
previously  given  for  the  example  problem  discussed  in  Section  IV.  The  results 
for  this  test  case  are,  in  general,  analogous  to  the  results  for  the  previous 
example  problem.  However,  if  one  examines  the  E(NORS)c  column  in  Table  XIV, 
he  finds  (1)  that  the  values  for  algorithms  4 and  5 are  not  monotonically 
decreasing  (at  a constraint  of  $ 78  M) ; and  (2)  at  two  constraints  ($  80  M 
and  $ 89  M) , algorithm  3 has  a value  lower  than  the  values  of  algorithms  4 
and  5.  This  can  be  explained  by  the  heuristic  nature  of  the  algorithms,  and 
by  the  fact  that  allocations  are  being  made  from  the  non-convex  regions  of  the 
item  functions,  (?,•„}  . Despite  an  occasional  small  reversal,  if  one  applies 
a little  decision  theory  (either  maximin  or  expected  value),  he  would  prefer 
to  use  algorithms  4 and  5 instead  of  algorithms  1,  2,  or  3 to  estimate  the 
minimum  value  for  E(N0RS)c. 

SAMPLE  OUTPUT  DISPLAYS 

Tables  XV,  XVI,  and  XVII  show  sample  printout  displays  for  the  test 
case.  Table  XV  shows  the  total  system  sum  of  expected  backorders  and  LCC 
corresponding  to  the  fifteen  Lagrange  multipliers. 

Table  XVI  is  a tabulation  of  the  KJ  variable  for  each  item  for  each 
of  the  15  Lagrange  multipliers.  The  definition  of  KJ  was  previously  discussed 
in  Supplement  A,  Table  XVII  shows  the  item  and  depot(s)  summary  printout  for 
a LCC  constraint  of  $ 89  M. 


TABLE  Xm 

COMPARISON  OF  ALLOCATION  ALGORITHMS  FOR  MINIMIZING  THE  SUM  OF 
EXPECTED  BACKORDERS  FOR  THE  TEST  CASE  WITH  SOME  u,  > 1 


SUM  OF 

EXPECTED  BACKORDERS 


1 2 


MAXIMUM  PERCENT  POSSIBLE 
REDUCTION  FOR  SOLUTION 


107  J 106.9 


76.8  76.1  76.1 


62.3  62.2 


52.4  52.4 


37.6  I 37.2  37.1 


26.6  26.1  25.1 


17.1  16.3  16.3 


10.2  10.1  10.1 


7.22  6.97 


6.07  I 4.71 


2.17  2.04  2.03 


0.4326 

0.0667 

0.0967 

0.9622 

M t ” rMl  TPI  Ml 

2.3776 

0.1913 

0.1096 

1.6663 

0.1660 

1.0499 

0.1697 

0.0221 

2.1996 

0.1623 

0.2271 

6.3138 

n::iE3i 

1.2068 

0.1706 

0.1706 

24.5728 

4.7079 

1.0027 

11.1667 

3.1642 

0.6919 

7.7146 

1.4632 

0.7391 

104.2 

103.8 

103.8 

78.2 

74.6 

74J 

62.6 

61.2 

61.2 

62.4 

61.7 

61.7 

37  J 

36.3 

26.4 

24.9 

24.9 

17.0 

16.2 

16.2 

10.2 

10.1 

10.1 

■o 

7.21 

MET1 

4.70 

4.69 

2.17  2.04  1 2.03 


"gfT—rr 


, 


j . 


TABLE  27 

TOTAL  SYSTEM  SUM  OF  EXPECTED  BACKORDERS  AND  LCC  FOR  THE 
TEST  CASE  WITH  SOME  U;  > 1 


K 

LA1BDAK) 

BT(K) 

LCC(K) 

1 

• 9152278E  +00 

. 24  09246E*0  3 

.70  59376EM8 

2 

• 1935  262E-04 

. 80  26935EV02 

•7175433E*98 

3 

.1258925E-04 

• 65  1 2454E»02 

• 7272958E+  9 8 

% 

• 1000800E-04 

• 53  81871E*0  2 

• 736 10  50  E+ 98 

5 

• 7 943282E-05 

* 48  25514E+  0 2 

• 74  495  9 IE*  9 8 

6 

• 6399573E-05 

• 33  0 2 227E*  02 

• 7657S91EO  8 

7 

.5011072E-05 

• 26 84659E* 0 2 

• 77 60575E+3  8 

• 

• 3931072E-05 

• 15  46232E+  0 2 

.80  18220E0  8 

9 

.3152278- -05 

• 1441885E*’  0 2 

• 80  46705E*3  8 

10 

• 2511886E-05 

• 10  6 3 240  E*- 0 2 

.81 72781E* 3 8 

11 

• 1935262E-05 

• 10  29U5E*02 

.8189371EH8 

12 

. 1554893E-85 

. 85  95054E*01 

• 8290134Ef08 

11 

• 1000000E-05 

* 5549581E*  0 1 

• 84  98272E06 

1% 

• 6309573E-06 

. 17  20900E*  01 

• 89458  31E>  3 8 

15 

• 3152278E-06 

• 12  24174E*’  01 

• 9066921Ef98 
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